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ABSTRACT 


A  time-varying  random  linear  system  model  is  hypothesized  for 
the  communication  medium.  The  fundamental  parameters  and  analysis 
techniques  for  the  linear  system  are  described,  and  applied  in  a  multiple 
alternative  communication  mode  with  Rayleigh  fading.  A  comparison  of 
several  receiver  types  is  then  effected.  Finally  some  measurement 
limitations  and  shortcomings  are  indicated. 
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1. 


Introduction 


The  detection  and  correct  identification  of  weak  signals  in  the  presence 
of  strong  noise  backgrounds!  intentional  and  otherwise!  is  an  important  pro¬ 
blem  in  many  fields.  In  particular,  it  has  held  the  rapt  attention  of  communica¬ 
tions  engineers  for  some  time  now.  Many  important  situations  and  problems 
have  been  analyzed  and  solved,  and  used  as  guides  for  dealing  with  other  more 
complicated,  but  similar,  cases,  and  much  of  this  wealth  of  information  and 
analysis  can  be  applied  to  the  problems  of  underwater  detection,  classification, 
and  communication. 

In  order  to  determine  the  applicability  of  existing  techniques  and 
solutions  to  improvement  of  ASW  system  performance,  as  well  as  to  assess 
the  needs  for  new  theory  and  techniques,  it  is  of  first  importance  to  know 
whether  it  is  valid  to  consider  that  source,  medium  (channel),  and  receiver 
can  be  distinguished  and  separately  characterized  in  a  one-way  system,  as  it 
can  in  the  usual  communications  application.  This  separation  is  desired  in 
such  a  way  that  the  receiver  output  due  to  an  arbitrary  source  signal  can  be 
predicted  from  knowledge  of  its  output  due  to  a  previously  applied  source  signal 
of  different  characteristics.  In  other  words,  can  a  model  of  the  system  be 
derived  which  will  permit  source  and  channel  to  be  independently  varied  in 
analysis,  or  are  the  two  so  inextricably  tied  together,  that  results  for  one 
signal  and  set  of  environmental  conditions  are  useless  for  predicting  results 
for  other  signals  and  other  channel  conditions?  We  doubt  that  the  latter  situation 
is  the  one  given  to  us  by  nature;  in  fact,  this  is  our  basic  working  assumption  in 
the  following. 

Even  when  separability  is  possible,  there  are  still  substantial  questions 
to  be  answered  concerning  the  channel;  e.  g. ,  whether  it  is  stationary,  deter¬ 
ministic,  etc;  what  is  the  detailed  structure  and  effect  on  transmitted  signals; 


and  last,  but  far  from  least,  what  measurements  should  be  made  in  order  to 
answer  these  questions?  In  the  following  we  shall  discuss  a  general  linear 
model  which  seems  suitable  for  a  first  try  at  characterizing  the  channel  and 
answering  these  questions.  There  are  several  reasons  for  choosing  a  linear 
model  of  the  channel.  First,  and  perhaps  most  important,  the  actual  medium 
is  believed  to  behave  linearly  for  limited  excitations.  Second,  the  analysis  of 
a  linear  model,  although  not  simple,  is  fairly  tractable  for  certain  quantities 
of  interest.  And  last,  optimization  of  the  performance  of  the  linear  model  is 
often  possible;  this  is  not  always  so  for  other  models  which  include  nonlinearities. 
The  importance  of  the  last  reason  can  not  be  overstressed,  for  it  often  provides 
methods  and  suggests  guidelines  to  be  considered  in  the  physical  application  to 
improve  performance. 

In  this  report  we  shall  present  some  of  the  fundamental  techniques 
available  for  studying  and  analyzing  random  time-varying  linear  systems  with 
both  qualitative  and  quantitative  results,  and  also  indicate  some  important 
limitations  on  measurements.  It  is  felt  that  this  work  will  furnish  an  adequate 
background  with  which  to  delve  deeper  into  additional  topics  of  particular 
interest. 


2. 


Linear  Model;  Generalitiea 


One  important  theme  running  through  virtually  all  current  communica¬ 
tion  work  ia  the  recognition  and  acceptance  of  the  transmission  medium  as  being 
a  linear  system,  although  perhaps  random  and  time-varying.  Since  the  trans¬ 
mission  medium  has  been  space  itself,  it  is  easily  observed  and  measured,  and 
tested  for  specific  applications.  Even  when  communication  is  over  long  ranges, 
through  ionospheric  disturbances,  or  via  moving  reflectors,  the  linear  model 
of  the  transmission  medium  is  an  excellent  one. 

Furthermore,  the  degree  of  randomness  of  the  linear  medium  is  usually 
relatively  small  in  comparison  with  the  deterministic  component.  That  is,  in 
most  cases  of  interest,  the  impulse  response  of  the  medium  h(t) (taken  time 
invariant  for  the  moment),  expressed  as  a  sum  of  a  deterministic  and  a  random 
part, 

h(t)  =  d(t)  +  r(t), 

has  negligible  energy  in  the  random  part  compared  with  that  in  the  deterministic 
part.  Mathematically  this  is  expressed  by 

^r^(t)  dt  <<  ^  d^(t)  dt, 

with  probability  near  unity.  Even  in  those  cases  where  the  impulse  response 
varies  with  time,  as  for  a  moving  transmitter  or  receiver,  or  for  communica¬ 
tion  through  the  ionosphere,  the  deterministic  part  of  the  impulse  response 
(that  part  which  can  be  predicted  from  recent  measurements)  dominates  the 
random  part. 


Integrals  without  limits  are  taken  over  the  range  of  non-zero  integrand. 
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The  communications  engineer  makes  great  use  of  the  deterministic 
components  of  both  radiated  signal  and  medium  impulse  response  by  utilizing 
matched  filters  (time  varying  in  some  cases)  at  the  receiver.  These  filters 
add  the  signal  portion  of  the  received  waveform  in  direct  proportion  to  the 
processing  time,  while  building  up  the  noise  only  as  the  square  root  of  the 
processing  time.  Hence,  for  long  processing  times,  the  output  signal-to-noise 
ratio  of  the  matched  filter  can  be  made  quite  large,  provided  that  the  receiving 
filter  can  be  kept  matched  all  this  time.  Mathematically,  if  s (t)  is  the  received 
signal  in  the  absence  of  additive  noise  (known  quite  accurately  to  the  receiver), 
and  n(t)  is  the  additive  (white)  noise  from  other  sources,  such  as  receiver 
noise,  environmental  noise,  or  jamming,  the  output  of  the  optimum  receiver  is 

r(T)  =  ^  s (t)  [  s(t)  +  n(t)]  dt, 

T 

where  T  is  the  processing  time  or  duration  of  the  signal.  This  quantity  has 
a  mean  value  of  E,  and  a  variance  N^E,  giving  a  (power)  signal-to-noise 
ratio  of 


2TW  | 


E  is  the  received  signal  energy,  is  the  double-sided  spectrum  level  of 
the  noise,  W  is  the  bandwidth  of  the  signal,  S  is  the  average  received  signal 
power  over  T,  and  N  is  the  average  received  noise  power  in  the  band  W.  The 
possibilities  of  improved  performance  are  well  evident  from  this  formula;  a 
large  TW  product  signal  can  combat  noise  very  well. 
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The  situation  for  underwater  transmissions  is  nowhere  near  as 
amenable.  The  validity  and  utility  of  the  linear  model  of  the  transmission 
medium  has  apparently  not  been  verified  to  any  high  degree  of  certainty. 

And  the  percentage  of  useful  deterministic  component  (much  less  than  the 
corresponding  figure  for  space  communication),  may  even  be  less  than  the 
percentage  of  random  component. 

To  further  complicate  the  situation,  the  radiated  "signal"  does  not 
maintain  a  high  degree  of  consistency  from  trial  to  trial,  whether  the  source 
be  friend  or  foe.  Even  for  a  friendly  submarine  repeating  a  run  at  the  same 
range,  speed,  bearing,  etc. ,  the  movement  of  water  over  the  transducers 
and  hull,  and  the  vibrations  of  the  auxiliary  machinery  and  the  submarine 
itself,  add  substantial  random  effects  to  the  intentionally  radiated  communica¬ 
tion  signal.  The  percentage  of  total  radiated  signal  which  is  deterministic  is 
apparently  not  known  with  any  great  degree  of  accuracy. 

In  order  to  be  able  to  combat  the  underwater  environment,  and 

corroborate  or  disprove  the  linear  model  of  the  medium,  some  quantitative 

measurements  should  be  undertaken.  The  guide  as  to  which  measurements 

may  be  furnished  by  considering  the  model  of  the  complete  communication 

system  given  in  Figure  1  below.  This  model  is  felt  to  be  a  fairly  reasonable 

characterization  of  the  actual  situation.  s(t)  is  that  portion  of  the  transmitted 

n(t) 


a(t)+a(t) 

h(T,t)  =  d(r,  t)  +  r(T,  t) 

y(t) 

•  w(t) 

1  ~ 

1 

Z 

i  ■ 

i 
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Linear  Time- 
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i 
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l 

p  meaium  - ^ 

Figure  1.  Model  of  the  Communication  System 
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waveform  which  ia  identical  from  trial  to  trial,  under  precisely  the  same 

environmental  conditions,  and  is  called  the  transmitted  signal,  a(t)  is  the 

* 

random  component  of  the  transmitted  waveform  and  is  unpredictable. 

The  medium  is  characterized  by  a  linear  time  varying  filter  with 
impulse  response  h(T,  t).  That  is,  h(T,  t)  is  the  response  (with  external  noise 
n(t)  *  0)  at  time  t  due  to  an  impulse  excitation  time  t  ago.  The  response  is 
measured  at  the  receiving  point  of  interest.  (If  the  medium  is  time  invariant, 
h(T,  t)  *  h(T). )  d(r,  t)  represents  the  deterministic  portion  of  the  impulse 
response  (predictable  from  recent  past  measurements)  and  r(T,  t)  the  random 
component.  The  output  y(t)  of  the  filter  (in  the  absence  of  external  noise) 
for  an  input  x(t)  is  given  by 


y(t)  = 


^  x(t-r)  h(r, 
j  x(t-T)  d(r. 


t)  dT 


t)  dT  + 


r(T,  t)  dT, 


and  is  composed,  likewise,  of  a  deterministic  and  a  random  part  (if  x(t)  is 

known).  Notice  that  the  random  component,  ^  x(t-T)  t(t,  t)  dT,  consists  of  a 
linear  operation  on  the  input  waveform  x(t),  and  represents  therefore  not 
merely  an  additive  noise  component,  but  more  like  an  averaged  multiplicative 
noise  effect  on  the  input  x(t).  If  x(t)  were  zero,  this  random  term  disappears--- 
quite  different  from  the  case  where  an  additive  noise  effect  is  still  present  when 
the  signal  is  absent. 


This  portion  of  the  waveform  is  not  useless.  Random  signals  for 
communication  and  detection  have  recently  attracted  the  attention  of  several 
investigators  [  l],  [2]  and  have  possibilities.  In  fact,  this  is  the  only  wave¬ 
form  on  which  detection  and  classification  can  be  attempted  when  the  submarine 
is  not  transmitting  a  known  communication  signal,  as  for  example,  while 
hovering.  However,  for  the  present,  in  order  to  retain  simplicity,  we  shall 
attempt  to  analyze  only  the  non- random  component  s(t),  and  save  the  generali¬ 
zations  for  the  future. 
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Since  the  ocean  is  not  ah  all-pass  linear-phase  network,  the  received 
signal  for  a  single  path  will  not  be  an  impulse  at  the  appropriate  delay  even 
if  the  excitation  were  an  impulse.  Rather  it  will  be  a  smeared  out  pulse 
centered  at  the  appropriate  range  delay.  For  example,  if  the  (time  invariant) 
linear  filter  has  a  frequency  domain  representation  approximated  by 


H(f) 


-  af  -i2irbf 
e  , 


f  >  0, 


the  impulse  response  would  be  approximately 


* 


2a _ 

2  2  2  ’ 
z  +  (2w r  (t-b) 

which  is  peaked  at  b,  with  a  width  of  a/ir  at  the  half  amplitude  points.  (The 
bandwidth  of  the  filter  is  21n  2/a  at  the  half  amplitude  points. ) 

Depending  on  the  particular  geometry,  range,  water  temperature,  etc. , 
there  may  be  one  or  more  major  contributions  to  the  medium  impulse  response. 
Thus,  for  a  situation  where  there  is  one  bottom  bounce,  the  impulse  response 
on  one  trial  may  appear  as  in  Figure  2. 


Figure  2.  Possible  Impulse  Response 
-And  if  fhe  difference  in  any  two  paths  is  small  enough,  their  pulse  outputs  will 
overlap.  But  there  is  still  no  need  for  an  excitation  which  is  narrower  than 
the  inverse  bandwidth  of  the  medium,  because  the  two  pulse  outputs  can  never 
be  resolved  any  finer  than  this  quantity  in  any  event. 

*Extremely  short  duration  excitations  are  not  necessary  to  measure  the 
impulse  response;  rather,  any  excitation  with  an  effective  duration  less  than 
approximately  a/w  in  the  example  above  is  acceptable  for  measuring  the 
impulse  response. 
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An  equivalent  way  of  characterizing  the  medium,  if  linear,  is  to 
measure  the  amplitude  and  phase  characteristics  in  response  to  pure  single 
tone  excitations.  The  function  givingths  information  is  called  the  instantaneous 
transfer  function  and  is  defined  by 

H{q ,  t)  =  j  h(r,t)  exp  (-i2irivr)  dr 
*  D(q,  t)  +  R(t)  ,  t) 

and  contains  both  deterministic  and  random  components.  If  a  sinewave 
A  sin  2irfjt 

$ 

is  transmitted,  the  received  waveform  (without  noise)  would  be  approximately 
A  |  D(fj,  t)  |  sin  [ 2irfjt  +  t)], 

provided  D(f^.t)  changes  very  slowly  over  times  approximately  l/ f ^  seconds 
in  duration. 

The  final  response  of  the  system  is  given  by 

[  s(t-T)  +  a(t-r)  ]  [  d(r,  t)  +  r(r,  t)  ]  dr  +  n(t) 

s(t-r)  d(T,  t)  dr 
s(t-r)  r(r,  t)  dr 
a(t-T)  d(r,  t)  dT 
a(t-r)  r(r,  t)  dr  +  n(t) 

*  w j (t)  +  w2(t)  +  w3(t)  +  w4(t). 

* 

With  noise  present,  this  expression  represents  the  mean  of  the 
received  waveform. 


•5 

♦J 

♦f 

*S 
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The  first  component  of  the  response.  Wj(t),  is  the  (desired)  well  be¬ 
haved  deterministic  component  of  the  received  waveform.. and  is  the  one  upon 
which  normal  reception  is  attempted. 

Component  w^  (t)  is  a  multiplicative  noise  term  resulting  from  the 
passage  of  a  deterministic  signal  through  a  random  filter.  Component  w^(t) 
is  deterministically  filtered  random  noise.  These  two  noise  components  may 
be  used  in  some  cases  to  improve  performance  above  that  of  w^(t)  alone. 
Whether  or  not  this  will  be  fruitful  depends  on  the  relative  magnitudes, 
compared  with  the  first  component. 

Component  w^(t)  is  noise  which  is  useless  and  should  be  suppressed 
as  much  as  possible.  Whether  or  not  this  goal  can  be  adequately  reached  again 
depends  on  the  relative  magnitudes  involved. 

It  is  impossible  to  measure  separately  any  of  the  components  {wk<t)} 
or  their  averages.  Rather,  they  come  in  an  inseparable  bundle,  and  statements 
about  each  type  can  be  obtained  only  through  measurements  of  s,  a,  d  and  r. 
Once  these  measurements  have  been  made,  however,  we  can  obtain  many 
statistical  quantities  of  interest  in  terms  of  them.  For  example. 


wk(t)  =  0,  k  =  2,  3,  4. 

w^(t)  =  sft-Tj)  s(t-T2)  r(Tlt  t)  r(T2,t)  dTj  dr2 

If  r(t)  is  stationary  such  that 


r(Tj,  t)  r(r2,  t)  =  R(t-Tj)  R(t-r2), 


we  have 


(t) 


[I 

[I 


s(t)  R(t)  dt  j 
G(f)  S*(f)  df 
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where  S(f)  and  G(f)  are  the  Fourier  transforms  of  s(t)  and  R(t)  respectively. 


Also 


W3  (t) 


a(t-Tj)  a(t-r2)  d^.tld^,  t)  dTj  dT2 


If  a(t)  is  stationary, 


aft-Tj)  a(t-r2)  =  A^-  rj, 

and 

w23(t)  =  Mr^  t2)  d(T  1#  t)  d(r2>t)  dTi  dT2 

=  j  <r  (f)  I  D(f,  t)  |  2  df 

where  or  (f)  is  the  Fourier  transform  of  A(t). 

Lastly,  using  the  independence  of  the  three  noise  terms  r(r,  t),  a(t), 
and  n(t),  we  have 

w2(t)  =  aft-Tj)  a(t-r2)  r(Tj,  t)  r(r2>  t)  dTjdT,,  +  n2(t) 

A^Tl"  T2^  R*t"1V  R^t-T2^  dTldT2  +  h2 

=  j  a  (f)  |  G(f)  |  2  df +  n2 
in  the  case  of  stationary  statistics. 


Combining  all  the  results  above,  we  have 


s(t-r)  d(r,  t)  dT 


2 

O2  [w(t)]  *  G(f)  S*(f)  df  +  J  a( f)  |  D(f,  t)  | 2  df 


£  »<f)  |  G(f)  I  2  df  +n2 


from  which  the  receiver  input  signal -to -noise  ratio  can  be  determined.  The 
major  contributions  to  the  receiver  input  noise  can  be  readily  evaluated  from 
this  equation.  The  areas  of  investigation  and  possible  improvement  are  then 
more  easily  ascertained.  Other  moments  of  interest  can  also  be  derived, 
provided  the  higher  order  correlation  functions  of  the  random  components 
are  known. 

The  statistics  described  above  are  function  of  time  in  general,  due 
to  the  non- stationary  character  of  the  random  fluctuations.  However,  it  is 
possible  to  perform  time  averages  on  these  statistics  and  relate  their  long 
term  averages.  Depending  on  the  duration  of  the  tests,  this  procedure,  a 
quasi- stationary  one,  may  be  the  only  alternative  open  for  computation. 

A  more  complete  statistical  description  of  w(t)  is  afforded  by  its 
p.  d.  f. ' s  and  not  just  its  moments.  However,  the  relation  of  the  p.  d.  f.  1  s 
of  w(t)  to  those  of  a,  d,  etc.  is  extremely  difficult  except  in  special  cases. 
One  of  those  special  cases  is  when  the  random  component  r(T,  t)  of  h(T,  t)  is 
negligible  in  comparison  with  the  deterministic  component  d(r,  t).  Then,  if 
no  signal  is  transmitted 


a(t-r)  d(r,  t)  dT  +  n(t) 
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where  we  recall  that  a(t)  it  the  random  emiasion  of  the  aubmarine.  If 
furthermore,  d(r,  t)  ariaea  due  to  juat  two  patha,  auch  aa 

d(T,t)  =  Pj  tir-AJ  +  P2  5(t-A2), 
where  P^,  P2<  A^,  A2  may  be  varying  slowly  with  time,  then 

w(t)  =  Px  alt-Aj)  +  P2  a(t-A2)  +  n(t). 

The  firat  order  p.  d.  f.  of  w(t)  is  then  given  by 

p(w- 1)  *  Pa(xj.  t-Aj;  x2.  t-A2)  Pn(w-PlXl  -  P2x2.  t)  dxjdx2 

where  the  second-order  p.  d.  f.  of  a  and  the  first  order  p.  d.  f.  of  n  are 
necessary.  Therefore  if  the  first  order  p.  d.  f's  of  w  and  n  and  the  second 
order  p.  d.  f.  of  a  are  measured,  they  may  be  substituted  above  and  checked. 
For  second  order  statistics,  or  higher,  of  w,  the  analytical  approach  is 
fraught  with  integrals  of  still  higher  order  p.  d.  f.  's  of  a  and  n.  Furthermore, 
if  d(r,  t)  does  not  contain  just  impulses,  the  problem  is  almost  completely 
intractable  except  in  very  special  cases. 

In  the  above,  we  have  tried  to  briefly  indicate  some  of  the  statistical 
quantities  of  interest  and  importance,  and  how  they  may  be  interrelated. 

In  the  following  sections,  we  shall  delve  much  deeper  into  particular  aspects 
of  the  communication  model  and  obtain  quantitative  results. 


3. 


Fundamental  Parameter «  of  a  Time  Varying  Linear  System 

In  eection  2  we  defined  the  filter  iropulae  response  of  a  general  time- 
varying  (random)  linear  system  as  h(T,  t),  the  response  of  the  system  at  time 
t  to  a  unit  impulse  applied  t  ago.  Also,  the  instantaneous  transfer  function 
of  the  system  was  defined  as 


H(f|.t)  *  J  h(r,  t)  exp  (-i2irivr)  dT 

and  was  indicated  to  measure  the  response  of  a  system  to  pure  tone  excitations. 
It  will  prove  convenient  to  define  two  more  equivalent  system  functions 
(characterizations)  for  later  use.  They  are  the  conditioned  spreading  function 
of  the  filter 

A(t,  f)  =  ^  h(T,  t)  exp  (-i2irft)  dt 
and  the  bi-frequency  function 


(Tl»f)  =  J  H(r),t)  exp  (-i2irft)  dt 
=  ^  A(t,  f)  exp  (-i2irqr)  dr 
=  jJh(T.t)  exp  [-i2w(iiT  +  ft)  ]  dr  dt. 


Suppose  h(T,  t)  is  substantially  confined  within  a  rectangle  L  by  D 
in  the  t,  t  plane;  that  is,  h(r,  t)  is  approximately  zero  outside  the  rectangle. 
L  is  called  the  "memory  time"  or  delay  spread"  of  the  filter,  while  D  is  the 
"duration"  or  "existence  time"  of  the  filter.  For  an  RC  filter,  for  example, 
D  *  to  and  L  ="  3RC. 
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Since  H(tj  ,  t)  is  the  Fourier  transform  on  r  of  h(T,  t).  H(r| ,  t)  is 
confined  to  a  strip  of  width  D  on  the  t-axis>  just  as  h(r,  t)  is.  Furthermore, 
for  t  in  this  strip,  suppose  H(T).t)  is  substantially  confined  within  a  band  of 
width  F  in  1^.  F  is  called  the  "bandwidth"  of  the  filter.  Then  H(ij,t)  is 
significantly  non-zero  only  within  a  rectangle  F  by  D  in  the  rj ,  t  plane.  For 
the  RC  filter  mentioned  above,  F  =  3/RC. 

Also,  since  A(r,f)  is  the  Fourier  transform  in  t  of  h(T,t),  A(t,  f) 
is  confined  to  a  strip  of  width  L  on  the  r-axis,  just  as  h(r,  t)  is.  Further¬ 
more,  for  T  in  this  strip,  suppose  A{t,  f)  is  substantially  confined  within  a 
band  of  width  B  in  f.  B  is  called  the  "frequency  spread"  of  the  filter.  Then 
A(t,  f)  is  significantly  non-zero  only  within  a  rectangle  L  by  B  in  the  t,  f  plane. 
Again  for  the  RC  filter,  B  =  0  (no  spreading). 

Lastly,  since  (r| »  f)  is  the  Fourier  transform  in  t  of  H(tj,  t),  and 
the  Fourier  transform  in  t  of  A(r,  f),  ^  (f) ,  f)  is  confined  within  a  rectangle 
F  by  B  in  the  t) ,  f  plane. 

The  four  diagrams  in  Figure  3  depict  "slices"  of  the  functions  in 
the  various  domains.  A  fundamental  property  of  Fourier  transforms  has 
been  used  extensively  here;  namely  if 

G(y)  3  g(x)  e  l2irxy  dx 

and  g(x)  has  extent  X  on  the  x-axis,  then  G(y)  can  change  significantly  in  a 
distance  no  smaller  than  on  the  y-axis,  provided  g  does  not  have  large 
alternating  values;  i.  e. ,  we  do  not  consider  "Buperdirectivity". 

A  couple  of  comments  should  be  made  about  Figure  3.  If  the  filter 
is  physically  realizable,  h(T,  t)  =  0  for  t  <  0,  and  the  rectangle  will  not 
extend  to  the  left  of  the  t-axis.  Also,  for  a  time  invariant  network  (and  many 
others),  D  =  co;  for  this  case,  the  wiggles  in  f  of  A(r,  f)  and  tj,  f)  occur 
infinitely  fast.  Figure  3  is  drawn  for  the  most  general  case.  . 
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Figure  3.  FILTER  FUNCTIONS 
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By  way  of  illustration  of  the  parameters  defined  above,  consider  that 

a  sine  wave  of  frequency  is  fed  into  the  filter  described  by  Figure  3. 

From  3(a).  the  output  would  be  significantly  non-zero  only  for  approximately 

the  D  seconds  indicated  on  the  t- scale  (D  +  L  seconds,  more  accurately. ) 

From  3(b).  there  would  be  no  output  at  any  time  unless  riQ  lies  in  the  pass 

band  F  of  the  filter.  And  even  when  v|  lies  in  this  region,  the  amplitude  of 

o 

the  output  sine  wave  is  scaled  by  the  factor  tfy  and  therefore  varies 

with  the  input  frequency  And  from  3(c),  the  spectrum  of  the  output 
signal  would  not  be  an  impulse,  as  for  the  input,  but  rather,  spread  out  on 
the  f- scale  over  a  range  B  wide. 

To  summarize, 

D  =  existence  time 
L  =  delay  spread 
F  =  bandwidth 
B  =  frequency  spread 
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4. 


Narrow-Band  Linear  Systems 


In  Section  3,  general  definitions  of  linear  system  characterizations 
were  given.  However,  these  characterizations  are  not  convenient  for  using  in 
situations  where  the  systems  are  narrow-band,  i.  e. ,  pass  only  a  small  band  of 
frequencies  about  some  high  center  frequency.  (In  terms  of  Figure  3(b)  or  3(d), 
this  means  the  rectangle  would  be  displaced  far  to  the  right  on  the  *1-  scale. ) 

Such  a  situation  can  arise  naturally  when  a  narrow  band  excitation  passes  through 
a  broad-band  medium,  for  reception  would  then  reasonably  consist  of  filtering 
out  only  that  portion  of  the  spectrum  where  the  signal  is  expected  to  lie.  Alter¬ 
nately,  a  piece  of  reception  equipment  may  be  inherently  narrow-band  and  able 
only  to  accept  a  restricted  region  of  input  signal  spectrum,  even  though  the 
input  is  broad-band.  In  this  section,  therefore,  we  specialize  the  results  of 
Section  3  to  narrow-band  systems,  and  indicate  some  of  the  more  usual  types. 

This  material  is  used  extensively  in  Section  5  for  deriving  performance  in  the 
presence  of  noise. 

We  had* 

H(f,  t)  »  ^  h(r,  t)  exp(-i2irfT)  dr 

as  the  instantaneous  transfer  function  of  a  linear  system.  For  an  input  i(t), 
the  output  r(t)  may  be  obtained  as 

r(t)  »  J  h(r,  t)  i(t  -  r)  dT 
*  ^  H(f,  t)  1(f)  exp(i2irft)  df 

where  1(f)  is  the  voltage  density  spectrum  (Fourier  transform)  of  the  input.  The 
latter  expression  in  the  frequency  domain  is  a  neat  generalization  of  the  usual 
time-invariant  result. 

$ 

We  use  the  more  conventional  symbol  f  in  this  section,  rather  than  q,  because 
there  is  no  chance  of  confusion  as  there  was  in  Section  3,  where  a  bi-frequency 
function  was  defined,  f  is  only  a  labeling  index,  and  is  not  a  fundamental  parameter. 
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Since  h(r,  t)  and  i(t)  are  real,  we  can  easily  show  that 


r(t)  *  Re  j  C  H(f,  t)  1(f)  exp(i2irft) 

(o 

(We  drop  all  irrelevant  scale  factors  in  order  to  facilitate  computations.  Our 
result  will  be  proportional  to  the  actual  output,  and  will  be  sufficient.)  Defining 


¥(f)  »  1(f)  [  1  +  sgn(f)] 

H+(f,  t)  =  H(f,  t)  [  1  +  sgn(f)] 


we  find 


r(t)  =  Re 


H+(f,  t)  'i'(f)  exp(i2irft)  df 


} 


(1) 


Now  for  a  narrow-band  system,  H+(f,  t)  will  be  significantly  non-zero  only  near 
one  frequency  f^  (for  all  t),  and  we  express  it  as 


H  (f,  t)  a  H  (f  -  f.  ,  t) 

+  LF  h 

where  HLF(f»  0  has  a  filter  characteristic  centered  about  zero  frequency- -a  low- 

)|( 

pass  filter.  Then  letting 


*(f)  =  Z(f  -  fh), 

we  find 

r(t)  =  Re  { exp(i2irf^t)  ^  t)  Z(f)  exp(i2rrft)  df} 

=  Re  {  expUZirf^t)  J  hLp(r,  t)  z(t  -  r)  dr} 

where 

hLF(T,  t)  =  J  Hj^p(f,  t)  exp(i2rrfr)  df 


The  "center  frequency"  of  ^(f)  need  not  be  f^.  Thus  Z(f)  need  not  be 

centered  at  zero  frequency.  In  fact,  Z(f)  -  Vr  (f  -  f  +  f  ),  where  f  is  the 

Lif  8  n  8 

input  signal  center  frequency. 
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and 


z(t)  ■  ^  Z(f)  exp(i2irft)  d£ 

»  $(t)  exp[i2ir(fg  -  f^t] 

where  $(t)  ia  the  "complex  low  frequency  envelope"  of  the  input  (Fourier 
transform  of  Therefore  the  complex  low  frequency  envelope  of  the 

response  r(t)  is  proportional  to 

E(t)  =  J  hLF^T’  ^  z ^  ’  7 *  dT‘  (2) 

h^p(r,  t)  is  the  complex  low  frequency  impulse  response.  A  block  diagram 
relating  the  complex  envelopes  is  shown  in  Figure  4. 


Linear 

Time- Varying 


z(t)  . 

hT  V(T.  t) 

E(t). 

7 

Jbr 

Figure  4.  Equivalent  Complex  Envelope  System 


Summarizing,  we  have 
r(t)  *  Re  {E(t)  exp(i2irf^t)} 

E(t)  =  J  1»lf(t,  t)  z  (t  -  t)  dr 
z(t)  »  C(t)  exp[  i2w(fg  -  £h)  t] 

i(t)  =  Re  {£  (t)  exp(i2irf  t)}  (3) 

s 


The  only  quantity  not  completely  tied  down  here  is  hT  _(t,  t).  We  now  proceed 

Lx 

to  relate  h^F  to  h:  let  i(t)  =  6(t  -  tQ);  then,  using  (1), 


r(t)  =  h(t  -  tQ,  t) 


»  Re 


t)  exp[  i2irf(t  -  t  )] 
o 


* 


See  [2],  Ch.  2, 
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*  Re  |  ^  HL,p(f  “  fh*  ^  '  *<,)]  df  j 

■  Re  {  exp[  ttirf^t  -  tQ)  J  HLF(t),  t)  exp[i2ini(t  -  tQ)]  dn  } 

■  Re  {  exp[  i2iri^(t  -  tQ) ]  hLF(t-t0'  *)) 

Letting  t-tQ  *  T  (to  eliminate  tQ),  we  obtain  finally 

h(r,  t)  =  Re  {hLF<T,  t)  expt^ir^T)}  (4) 

(which  checks  with  intuition).  This  relation  allows  us  to  obtain  b^F  easily  from 
h.  (It  is  interesting  to  note  that  we  cannot  derive  (4)  directly  from  (2)  because 
f#  ■  <o  for  the  impulse,  and  we  get  an  uninterpretable  result. ) 

Special  Cases 

In  (2),  h  (t,  t)  is  a  complex  function,  and  furthermore,  can  be  random 
Lx 

or  deterministic.  For  example  if 
hLF(T'  t]  =  G(t)  6(T) 

where  G(t)  is  a  complex  random  process,  then 
HLF(f.  t)  =  G(t) 

which  is  independent  of  fi  This  describes  a  frequency  non- selective  medium  with 
random  gain  G(t).  The  complex  envelope  of  the  received  signal  is,  from  (2), 

E(t)  =  G(t)  z(t). 

The  block  diagram  of  this  system  (a  special  case  of  Figure  4)  is  depicted  in 
Figure  5. 

G(t) 

iii - ^ 

Figure  5.  Frequency  Non-Selective  System 
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From  (4), 

h(T,  t)  *  6(t)  Re  { G(t)} 


and  the  impulse  response  of  the  system  is  merely  a  time  varying  impulse  with 
no  delay--a  variable  gain  system.  For  a  time-invariant  system,  G(t)  *  V  and 
the  usual  simpler  results  follow. 

An  alternate  form  of  (2)  is  given  by 
E(t)  =  J  HLF(f,  t)  Z(f)  exp(i2irft)  df 

Now  if,  for  all  t,  H^p(f,  t)  is  roughly  independent  of  f  in  the  range  where  Z(f) 
is  non-zero,  then 

E(t)  =  HLF(fs  -  V  t}  Z(t) 

and  we  have  frequency  non- selective  fading.  Thus  we  obtain  the  reasonable 
result  that  H^p(f,  t)  need  not  be  independent  of  f  for  all  f,  but  only  over  the 
signal  bandwidth,  in  order  to  realize  frequency  non- selective  fading. 

As  a  second  example,  suppose 
hLF(T,  t)  =  G(t)  a(T). 

Then 

HLF(f,  t)  *  G(t)  0(f) 

where  &({)  is  the  Fourier  transform  of  a(r).  This  is  a  frequency  selective 
(random)  system.  There  follows  that  the  complex  output  envelope  is  given  by 

E(t)  =  G(t)  J  a(r)  z(t  -  r)  dT 

=  G(t)  ^  Q.(f)  Z(f)  exp(i2irft)  df 

The  block  diagram  of  Figure  6  depicts  this  system,  in  terms  of  the  complex 
envelopes. 
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Linear 

Time-Invariant 


z(t) 


*  a(r) 


G(t) 

I 

-s®- 


E(t) 


Figure  6.  Frequency  Selective  System 

The  input  signal  passes  through  a  filter  of  characteristic  #(f)  fvhich  may  be 
random)  followed  by  a  "gain  control"  G(t).  For  G(t)?l,  we  have  the  time-dispersive 
medium.  For  more  general  G(t),  we  have  a  time-  and  frequency-dispersive 
medium  (L,  B  >  0).  If  a(T)  is  not  random,  an  important  statistical  quantity  of 
h^p(T,  t)  is  the  correlation  function  of  the  medium  fading  G(t): 

G^)  G*(t2)  ®  d>(t1  -  t2) 

for  stationary  fading. 


As  a  last  example,  consider 

hLF(T’  Z  Gk<1:)  ak(T) 
k 

for  which 

hlf«  *»  *  l  Gk<*> 

k 

and 

E(t)  *  ^  Gk(t)  J  Z(f)  exp(iZirft)  df 

k 

which  is  a  parallel  combination  of  the  building  block  depicted  in  Figure  6.  For 
ak(T)  =  6(t  -  rk), 

E(t)  *  ^  Gk(t)  z(t  -  Tfc) 
k 

which  is  the  multipath  spreading  model  often  met  in  practice. 

It  is  thus  seen  that  Figure  4  depicts  a  very  general  situation  for  narrow- 
band  communication;  we  shall  analyze  a  special  case  of  it  in  the  next  section. 
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5. 


Analysis  of  a  Fading  Medium 

In  this  section,  results  of  an  early  study  of  the  performance  of  several 
different  receiving  structures  are  summarized.  The  pertinent  derivations  are 
presented  in  Appendices  A,  B,  and  C. 

The  model  to  be  studied  is  still  the  one  depicted  in  Figure  1  of  Section  2. 
However,  several  specializations  have  been  made.  First,  there  is  no  random 
component  in  the  transmitted  waveform:  a(t)  =  0.  Furthermore,  s(t)  should  more 
properly  be  replaced  by  s^(t)  where  1  <  k  <  M;  that  is,  the  transmitted  signal  is 
one  of  M  (equi-probable)  alternatives.  Second,  h(T,  t)  will  be  of  the  form  des¬ 
cribed  in  Section  4,  i.  e. ,  narrow-band.  Specifically,  the  first  special  case 
described  there,  Figure  5,  will  be  utilized,  and  G(t)  will  be  assumed  to  be  a 
complex  Gaussian  process.  Thus  the  fading  is  assumed  to  be  Rayleigh  distri¬ 
buted  and  non-frequency  selective.  Third,  the  additive  noise  n(t)  is  assumed  to 
be  white  Gaussian.  These  restrictions  have  been  introduced  mainly  to  simplify 
the  analysis;  however,  they  do  represent  a  very  reasonable  state  of  affairs. 

The  three  receiving  systems  to  be  compared  here  are  the  M-ary 

)|CJ|C 

Weighted  Radiometer,  the  so-called  M- M-ary  System,  and  the  M-ary  Linear 

Filtering  System.  The  first  system  is  the  optimum  system  (from  the  standpoint 

of  minimizing  bit  error  probability  using  an  M-ary  alphabet)  when  the  signal 

energy-to-noise  power  density  ratio  is  small  over  a  coherence  interval,  i.  e.  , 

S  T 
A  o 

when  — ^ —  <<  1.  r  is  the  length  of  time  over  which  it  is  possible  to  reasonably 
o 

predict  the  amplitude  and  phase  of  the  medium  output  given  the  input.  In  terms 
of  the  correlation  function  of  the  medium  fading  introduced  in  Section  4,  we  define 
t  ,  the  coherence  interval,  such  that 


*[3] 
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I  ♦(t  )  |  *  |  <Mo). 

The  M-M-ary  method  of  reception  is  related  to  the  Weighted  Radiometer, 
and  produces  only  slightly  higher  bit  error  probabilities.  The  third  technique  con¬ 
sists  of  the  optimum  system  for  use  in  a  non-fading  channel,  and  has  been  included 
in  this  early  investigation  to  illustrate  the  fact  that  it  is  extremely  hazardous  to 
presume  that  a  system  optimized  for  combatting  one  type  of  disturbance  in  the 
transmission  medium  will  provide  near -optimum,  or  even  good  performance, 
against  another  type  of  disturbance. 

It  has  been  demonstrated  for  a  non-fading  medium  that  as  M— *  oo,  for 
both  phase-coherent  and  phase-incoherent  receivers,  the  probability  of  error 
— ^0  if  the  source  information  rate  H  is  less  than  the  capacity  of  the 
infinite  bandwidth  medium: 

(S  S 

1  +  TTW  *  18  6  if  bits/sec- 

o  o 

where  N  is  the  single- sided  noise  power  density  level,  and  S  is  the  average 
o  & 

signal  power.  Since  such  systems  are  theoretically  almost  100  %>  efficient  in 
the  usage  of  the  channel,  at  least  in  the  non-fading  environment,  their  wide¬ 
spread  usage  for  other  environments  is  suggested.  However,  some  method  of 
combatting  or  tolerating  fading,  fast  or  slow,  must  be  incorporated,  for  present 
designs  perform  miserably  in  the  presence  of  fading,  as  will  be  shortly  demon¬ 
strated. 


Basically,  it  is  anticipated  that  M-ary  systems  designed  for  fading 
should  be  near  optimum  in  channel  usage.  This  feeling  is  based  on  the  fact 
that  adequate  communication  through  a  fading  channel  is  possible  only  if  some 
averaging  of  instantaneous  medium  fluctuations  is  effected,  thereby  deriving 


[3,5,  6,  7] 
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the  benefits  of  the  average  medium  characteristics.  M-ary  signals,  by  their 
very  nature,  tend  to  occupy  sizeable  portions  of  the  time  axis,  and  experience 
multiple  fades  during  transmission  of  a  single  message  of  lg  M  bits.  There¬ 
fore  we  have  concentrated  in  this  section  on  modifications  of  the  basic  coherent 
M-ary  receiver  in  an  attempt  to  combat  fading  and  achieve  reliable  efficient 
performance;  it  is  believed  this  has  been  accomplished. 


We  now  embark  on  a  detailed  description  of  the  M-ary  systems  con¬ 
sidered  in  this  report.  If  the  usual  M-ary  system  is  used  in  a  fading  environment, 
and  orthogonal  signals  are  used,  a  simple  generalization  of  some  past  work* 
shows  that  the  probability  of  correct  decision  per  (M-ary)  word  is  given  by 
M-  ! 


whe  re 


Pc  =  TT  [*  +  n(l  +  0)  ] 
n=i 

0  =  —  j  |l  -  l«Mu)  |  2  du 


and  <f  is  the  correlation  function  of  the  medium  fluctuations.  For  a  linear 
correlation  function 


<(>(u)  =  1 


2  r 


u  <  2t 


it  is  readily  demonstrated  that 


0  = 


S  T 
a 

N 


1-1  — 
6  r 

o 


T  <  2t 


2t  i  ,  r 

—2.  i  .1  _2. 

T  3  T 


T  >  2t 
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Thus  even  as  T. — ►»,  0  saturates  at  2  — —  .  and  the  error  probability  does 

o 

not  tend  to  zero.  This  is  the  severe  limitation  of  the  standard  M-ary  technique. 

An  M-ary  system  of  great  generality  and  potentiality,  not  suffering  the 
above  limitations  is  realized  by  transmitting  a  sequence  of  one  of  q(<  M)  wave¬ 
forms  n  times,  and  employing  linear  weighting  to  realize  M  different  possibilities. 
For  example,  with  q  =  2,  and  Hadamard  matrix  weighting  of  size  M,  M  orthogonal 
signals  are  realized;  yet  only  two  possible  signals  need  be  generated.  Another 
special  case  is  realized  when  q  =  M;  this  case  is  called  M-M-ary  transmission. 
Here,  one  of  M  waveforms  is  repeated  several  times,  and  the  outputs  of  corres¬ 
ponding  matched  filters  appropriately  combined  to  yield  the  set  of  M  decision 
variables.  A  timing  diagram  of  the  waveforms  in  such  a  system  is  indicated 
below  in  Figure  7. 


Figure  7.  Timing  Diagram 

t  is  the  pulse  width,  A  is  the  pulse  spacing,  and  T  is  the  message  duration. 
The  probability  of  correct  decision  (per  M-ary  word)  is  derived  in  Appendix  B 
and  is,  for  v -th  envelope  combining,  and  for  A  >>  r  ,  r  <<  r  , 
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where 


Pc  "  J  ♦(x)  [ax  +  0]  dx* 
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yv  = 


r  iI  +  i 


\/r  (v  +  i)  -  r2(|  +  i) 


4>(x)  =  (2ir)  exp  (-  xZ/2) 


«(x) 


J  <j»(y)  dy. 


-oo 


The  limitations  on  A  and  r  above  ensure  that  independent  fading  occur  for 
each  pulse,  and  negligible  fading  occur  within  a  pulse.  These  assumptions 
allow  averaging  of  medium  characteristics  over  many  coherence  intervals. 


For  a  peak  power  limitation  on  the  transmitter,  it  may  be  impossible 
to  realize  the  large  pulses  required  by  the  M-M-ary  system.  In  such  a  case, 
the  space  between  pulses  should  be  filled  with  signal  power,  and  efficient  use 
made  of  such  waveforms.  For  the  limiting  case  of  peak- equal-average  power, 
3^  *  S^,  a  continuous  transmission  of  signal  power  over  the  message  interval 
T  is  dictated.  If  also 


S  r 
a  o 

N 


<<  1, 


** 


the  optimum  receiver  is  known  ,  and  is  called  the  weighted  radiometer.  One 


ihis  integral  is  tabulated  for  a  few  selected  values  of  a  and  8  in 
TR-63-3-BF,  "A  Multiple  Alternative  Error  Integral",  by  Albert  H.  Nuttall, 
April  30,  1963. 

**[  4] 
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of  M  possible  signals  is  generated  at  the  transmitter;  at  the  receiver,  M 
branches,  of  which  the  j-th  is  depicted  below  in  Figure  8,  are  constructed 
and  the  M  decision  variables  {D^}  compared,  a  decision  being  made  in  favor 
of  the  largest  variable. 


Figure  8.  Weighted  Radiometer 

The  bandwidth  of  H(f)  is  chosen  so  as  to  match  the  spreading  of  the  medium. 
The  probability  of  correct  decision  is  derived  in  Appendix  A  and  is 

P  ^  <(>(x)  {  x  +  3)  dx 

where 


For  a  linear  correlation  function. 


o  o 


St  ST 

&  O  3, 

Now  since  — —  <<  1  by  assumption,  we  see  that  — —  >>  1  is  required  for 
o  o 

good  performance.  Thus  extremely  large  message  durations  are  required 

for  this  system  to  achieve  a  high  degree  of  reliability. 

To  illustrate  the  degree  to  which  these  M-ary  systems  (with  M  >  2) 
are  superior  to  binary  systems,  four  binary  systems  have  also  been  evaluated: 
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the  binary  Weighted  Radiometer,  the  2-2-ary  system,  a  binary  linear  filtering 
system,  and  a  delay  autocorrelation  technique.  The  first  three  of  these  systems 
are  derived  from  the  M-ary  systems  with  M  =  2.  The  weighted  radiometer 
system  is  the  optimum  binary  system  for  low  signal-to-noise  energy  density 
ratios.  The  fourth  binary  system  represents  one  of  the  many  different  ways 
that  fairly  good  performance  can  be  attained  with  somewhat  less  complex  equip¬ 
ment  than  is  required  for  an  optimum  system;  it  is  analyzed  in  Appendix  C. 

Because  of  the  widespread  reliance  on  binary  digital  representations 
of  information,  performance  is  measured  by  the  bit  error  probability,  P^,  for 
each  system,  i.  e.  ,  the  long  term  average  probability  of  a  single  binary  information 
digit  being  assigned  the  wrong  value. 


The  bit  error  probability  for  each  of  the  four  binary  systems  is  plotted 


in  Figure  9  as  a  function  of  the  quantity  p  = 


where  H  is  the  information 


o 

rate  in  bits  per  second.  These  curves  are  applicable  to  the  case  where  the 
short-term  (peak)  power  limitation,  S^,  is  equal  to  the  long-term  average 
power  capability  of  the  transmitter,  S^,  Also  it  has  been  assumed  that  the 
fading  is  such  that  t^H  *  for  each  system.  The  most  striking  characteristic 
of  Figure  9  is  the  extremely  poor  performance  of  the  linear  filtering  system; 
although  optimum  for  use  in  a  non-fading  medium  (under  phase  incoherent 
operating  conditions),  its  performance  in  the  presence  of  fading  is  quite  un¬ 
satisfactory.  Another  feature  of  interest  is  that  the  weighted  radiometer,  which 
is  the  optimum  binary  system  for  small  values  of  p,  is  the  best  system  of  the 
four  for  all  values  of  p. 


Turning  to  the  multiple  alternative  systems,  the  performance  of  the 
M-ary  Weighted  Radiometer,  M- M-ary  and  M-ary  Linear  Filtering  systems  is 
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indicated  by  the  curves  in  Figures  10,  11,  12,  and  13.  The  weighted  radio¬ 
meter  has  been  evaluated  for  any  value  of  M,  under  the  constraint  of  peak- 
equal-average  power.  Curves  are  drawn  in  Figure  10  relating  the  bit  error 
probability  to  p,  for  M  =  2,  4,  8,  and  16,  assuming  that  the  fading  is  such 
that  t^H  *  ■—  .  As  expected,  the  bit  error  probability  is  reduced  significantly 
as  the  value  of  M,  the  alphabet  size  is  increased.  As  M  tends  to  infinity,  if  the 
M-th  order  Gaussian  approximation  involved  in  the  evaluation  of  error  proba- 
bility  is  valid,  then  the  limiting  performance  of  this  system  is  indicated  by 
the  dashed  curve  in  Figure  10.  This  suggests  (but  doesn't  rigorously  prove) 
that  the  limiting  performance  of  an  uncoded  M-ary  weighted  radiometer  achieves 
approximately  13  percent  of  the  medium  capacity,  C^,  under  the 
constraint. 

Analogous  curves  are  plotted  in  Figure  11  for  the  M-M-ary  system, 

under  the  same  constraints.  These  curves  indicate  that  the  M-M-ary  system 

is  approximately  0.  6  db  poorer  than  the  M-ary  radiometer  for  any  value  of  M. 

The  curves  in  Figure  12  represent  the  performance  of  the  M-M-ary  system 

when  essentially  no  constraint  on  peak  power  exists,  i.  e. ,  S  >  10  S  .  Under 

p  —  a 

these  conditions  a  method  of  averaging  over  many  independent  fading  coherence 
intervals  can  be  implemented  with  the  M-M-ary  system  to  improve  its  perfor¬ 
mance  significantly,  as  indicated  by  the  curves.  Again,  if  the  M-th  order 
Gaussian  approximation  is  valid  in  the  limit,  the  limiting  performance  Of  the 
M-M-ary  system  (indicated  by  the  dashed  line)  achieves  approximately  56 
percent  of  the  medium  capacity,  C  .  Although  time  has  not  permitted  the 
evaluation  of  the  M-ary  weighted  radiometer  under  the  no  peak  power  con¬ 
straint,  its  slightly  better  performance  than  the  M-M-ary  system  under  the 

S  »  S  constraint  suggests  that  the  former  system  may  be  capable  of  achieving 

pa 

nearly  100  percent  of  the  medium  capacity,  even  in  the  presence  of  fading. 

”  ^Although  this  approximation  is  believed  to  be  valid,  time  has  not  permitted 
the  completion  of  a  rigorous  proof. 
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Figure  10.  Performance  of  M-ary  Weighted  Radiometer  With  a 
Peak -Equal-Average  Power  Limitation 
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Figure  12.  Performance  of  M-M-ary  System  With  No  Peak 

Power  Limitation 
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The  M-ary  Linear  Filtering  Performance  curves  are  plotted  in 

Figure  13  for  both  peak  power  constraints,  i.  e. ,  S  =  S  ,  and  S  >  10  S  . 

p  a  p-  a 

The  dashed  curves  correspond  to  the  constraint  S  ■  S  ,  and  the  solid  curves 

p  a  ! 

apply  when  S  >  10  S  .  As  before,  the  fading  is  characterized  by  t  H*  -r— . 
p“  4  o 

Although  increasing  the  alphabet  size  does  provide  some  improvement  in  the 
system,  basically  it  is  unable  to  cope  with  a  fading  medium. 

The  performance  of  each  of  these  three  M-ary  systems  is  plotted  in 
Figure  14  for  M  ■  16,  under  the  peak-equal-average  power  constraint,  and 
assuming  tqH  =  ~  . 

The  curves  presented  in  Figures  9-14  suggest  strongly  that  higher 
order  alphabets  hold  the  key  to  satisfactory  performance  at  high  information 
rates.  Further,  the  use  of  systems  optimized  for  a  non- fading  medium  can 
produce  extremely  poor  performance  in  a  fading  medium.  Therefore,  it  iB 
of  prime  importance  to  place  major  emphasis  on  the  determination  of  optimum 
system  structures  for  fading  media,  and  the  development  of  implementable 
block  diagrams  which  exploit  these  structures  with  a  large  alphabet. 


* 

These  are  based  on  the  work  reported  in  [8]. 
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Figure  13.  Performance  of  M-ary  Linear  Filtering  System,  With 
and  Without  a  Peak  Power  Limitation 
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Figure  14.  Performance  of  Weighted  Radiometer  M-M-ary,  and  Linear 
Filtering  Systems  With  M  =16,  S  =  S  ,  and  t  H  =  rg-  .  -37 


6. 


Some  Measurement  Considerations 


Many  of  the  results  presented  earlier  have  presumed  knowledge  of 
the  pertinent  medium  responses  and  average  behavior.  In  this  section,  we 
shall  consider  some  of  the  problems  in  attempting  to  determine  these  quantities 
from  measurements. 


6.  1  Direct  Measurement  of  a  Non- Random  Time- Variant  System 

Let  us  consider  the  direct  measurement  of  h(T,  t).  If  we  excite  this 

filter  with  an  impulse  at  time  t  ,  the  response  is  h(t  -  t  ,  t)  and  we  have  a 

o  o 

"slice"  of  h  along  a  45°  line  in  the  r,  t  plane.  Since  this  function  has  duration 

L  seconds  on  the  t-axis  (see  Figure  3(a)),  we  cannot  re-excite  the  filter  with 

an  impulse  to  get  another  "slice"  until  L  seconds  later  than  the  first  excitation; 

otherwise,  the  outputswould  overlap,  thereby  generating  ambiguities.  Thus  we 

can  measure,  in  succession,  h(t  -  t  ,  t),  h(t  -  t  -  L,  t),  h(t  -  t  -  2L,  t),  etc. 

o  o  o 

In  order  that  we  not  be  missing  any  significant  changes  in  h(r,  t)  by  this  technique, 

there  should  be  little  change  in  the  sequence  of  h-functions  so  obtained.  That  is, 

h(t  -  t  ,  t)  delayed  by  L  seconds  should  look  like  h(t  -  t  -  L,  t)  for  all  t,  t  where 
o  o  o 

the  h-functions  are  not  small.  Thus  we  want 


h(t  -  t  -  L,  t  -  L)  S'  h(t  -  t  -  L,  t) 
o  o 

for  all  t,  tQ  giving  h  arguments  in  the  L  by  D  rectangle  of  Figure  3(a).  Alter¬ 
nately  this  can  be  written  : 


h(T,  t  -  L)  =  h(T,  t) 


in  the  L  by  D  rectangle.  Since  for  fixed  r,  h(r,  t)  can  change  in  a  distance 
seconds  on  the  t-axis  (see  Figure  3(a)),  this  requirement  of  slow  change  is 
satisfied  only  if 


B 


or 


L  < 


1 

B 


BL  <  1. 
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When  this  condition  is  satisfiedi  interpolation  from  the  45°  { h(t  -  t  -  k  L,  t)} 
slices  to  the  entire  h(r,  t)  function  for  all  t,  t  is  then  possible. 

If  BL>  >  1,  it  is  still  possible  to  get  slices  of  h(r,  t)  as  above.  How¬ 
ever,  the  slices  may  differ  significantly,  and  it  is  no  longer  possible  to  inter¬ 
polate,  and  obtain  the  entire  function  h(r,  t)*.  This  is  the  limitation  of  over¬ 
spread  (BL  >  1)  filters. 

Notice,  however,  that  if  the  filter  could  be  "reset"  to  an  earlier  time 
(as  for  example  with  a  network  with  a  switch),  we  could  interleave  additional 
slices  of  the  h(T,  t)  function,  each  set  of  slices  necessarily  being  separated  by 
L  seconds  to  avoid  overlap,  of  course.  This  resetting  procedure  could  be  con¬ 
tinued  until  the  density  of  slices  was  sufficient  to  interpolate  to  the  entire  function. 
Thus  if,  for  example,  the  impulse  response  were  desired  known  over  the  inter¬ 
val  of  D  seconds  on  the  t-axis,  a  processing  time  of  BLD  seconds  would  be 
necessary  to  evaluate  the  entire  h(T,  t)  function. 

It  might  be  supposed  at  first  that  the  limitation  BL,  <  1  is  due  to  the 
particular  method  of  measurement  of  the  system,  namely  impulse  excitation. 
However,  we  run  up  to  the  same  identical  relation  when  we  consider  the  direct 
measurement  of  H(t),  t).  If  we  excite  this  filter  with  a  pure  tone  e^^o* 
starting  at  time  tQ,  the  output  is,  for  t  >  t^, 

t-t 

e1  ^o  J  h(r,  t)  e  0  dr 

-«o 
i2wT)  t 

=  e  °  H(r,  ,  t) 
o 


This  situation  is  similar  to  sampling  a  function  band-limited  to  W  cps, 
slower  than  the  rate  of  2W  per  second.  It  is  not  possible  to  interpolate  to  the 
entire  function  from  the  "sparse"  samples. 
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for  t  -  t  sufficiently  large,  i.  e. ,  outside  the  region  of  L  seconds  on  the  T-axis. 

Thus  by  exciting  the  filter  sufficiently  long  ago,  we  can  evaluate  H  on  any  line 

i]  ,  tin  the  t),  t  plane.  However,  in  order  to  be  able  to  evaluate  H  at  other 
o 

values  of  tj,  we  must  excite  the  filter  simultaneously  with  a  set  of  pure  tones. 
(This  procedure  is  the  exact  dual  of  the  measurement  of  h(T,  t).  There,  we 
excited  with  a  set  of  impulses  in  the  time  domain.  Here,  we  excite  with  a  set 
of  impulses  in  the  frequency  domain. )  In  order  not  to  miss  significant  changes 
of  H(t),  t)  in  t|,  we  must  space  the  pure  tones  by  no  more  than  7*  cps  (see 
Figure  3(b)).  If  the  pure  tones  are  at  frequencies  {q  ,  the  response,  after 
sufficient  time,  is  approximately 


I 


i2irq  t 

e  *  H(r,k,  t). 


In  order  to  be  able  to  evaluate  each  H(iik,  t)  separately,  significant  changes 
in  H(T)k>  t)  must  not  occur  in  less  than  L  seconds;  otherwise  it  would  be 
impossible  to  separate  the  individual  filter  output  toneB.  From  Figure  3(b), 
we  see  this  is  possible  only  if 

or 

BL  <  1. 


If  this  condition  is  satisfied,  H(t],  t)  may  be  interpolated  from  the  slices 
{H(V  t)}.  From  the  complete  H-function,  h(r,  t)  could  of  course  be 
obtained  by  Fourier  transformation. 


Both  measurement  methods  above  require  BL.  <  1  in  order  that 
complete  measurement  of  the  filter  be  possible.  Undoubtedly  if  other  test 
procedures  for  the  direct  evaluation  of  the  A-  and  ^.functions  were  proposed, 
they  too  would  require  BL  <  1  for  unambiguous  filter  evaluation.  The  re¬ 
striction  BL  <  1  is  a  fundamental  restriction  on  ability  to  measure  completely 
a  filter's  characteristics.  Furthermore,  we  repeat  that  this  filter  had  no 
random  variation  whatsoever;  random  variations  can  only  generate  tighter 
restrictions. 
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6.2 


Measurement  of  a  Random  System 


Measurements  of  the  two  components  d(T,  t)  and  r(T,  t)  of  a  general 
system  h(T,  t)  could  be  accomplished  (for  a  particular  geometry,  thermocline, 
known  range,  etc. )  by  using  an  explosive  source,  and  measuring  the  response 
at  the  desired  receiving  site,  provided  the  external  noise  n(t)  is  small  in 
comparison  with  the  random  component  r(T,  t).  The  signal  source  (in  the 
submarine)  must  be  turned  off  for  this  measurement.  In  fact,  the  source 
should  probably  be  removed  from  the  area  because  its  random  output  affects 
the  measurements. 

Since  an  explosive  source  is  itself  random,  the  response  of  the 
medium  will  contain  a  random  component  due  to  the  source.  However,  if 
A  is  the  effective  duration  of  the  "impulse",  estimates  of  d  and  r  up  to  frequencies 
less  than  l/A  will  be  okay.  Equivalently,  estimates  of  D  and  R  will  be  valid  for 
f<  1/A. 

These  measurements  should  be  repeated  as  rapidly  as  possible  in 
succession  a  number  of  times,  before  any  important  parameters  have  changed, 
such  as  temperature  of  the  water.  Of  course,  a  new  explosion  and  measurement 
must  not  be  undertaken  before  all  the  multiple  echoes  from  the  different  possible 
bottom  bounce  reflections  (multipaths)  have  been  received  from  the  past  explosion, 
as  discussed  in  Section  6.  1.  These  measurements  should  be  made  at  a  great 
enough  distance  from  the  source  of  the  explosion  in  order  to  avoid  near  field 
effects,  and  yet  not  so  far  away  that  the  external  noise  n(t)  dominates  r(t).  The 
tests  should  be  monitored  by  measuring  the  intensity  of  the  explosion  as  close 
to  the  source  as  reasonable  and  practical. 

In  order  to  ascertain  the  degree  of  non-stationarity  of  h(r,  t)  (more 
particularly  d(r,  t)),  the  series  of  measurements  of  h(T,  t)  should  be  divided 
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into  groups,  the  members  of  each  of  which  have  been  determined  under  sub¬ 
stantially  the  same  conditions.  The  size  of  each  group  should  be  large  enough 
so  as  to  be  able  to  determine  d(T,  t)  accurately  (by  averaging  the  measure¬ 
ments))  and  yet  not  so  large  that  time  varying  items(  such  as  water  temperature, 
have  made  any  significant  change  during  that  time.  For  example,  suppose  time 

is  separated  into  intervals  {T.}  such  that  during  each  interval  T.f  h(T,  t)  S'  h.(r) 

J  J  J 

for  t  *  T..  (Also,  of  course,  d(r,  t)  S'  d.(r)  for  t  «  T..  This  is  the  mathematical 
J  J  J 

way  of  saying  that  the  time  variant  parameters  change  little  during  T.. )  Then 

J 

suppose  a  series  of  measurements  {  w.,  (t)}  of  the  impulse  response  were  re- 

JK 

corded,  where,  for  j  fixed,  t.,  «  T.,  k  =  1,  2,...,  n.  That  is,  n  measurements 

jk  j 

of  the  impulse  response  were  made  during  the  time  interval  T^.  Then  an  approx¬ 
imation  to  the  deterministic  component  of  the  impulse  response  is  realized 
according  to 

n 
1 

d. 

n 

k=l 


VT)?  :  2  VT+V- 


The  shift  t.^  aligns  all  the  impulse  responses  at  "zero"  time,  which  is  chosen 
for  convenience.  As  j  progresses  through  its  values,  local  averages  of  d(r,  t) 
in  intervals  of  length  T^  on  the  t-scale  are  obtained.  That  is,  the  quantity 


Tj  ^ 

J  rp 


d(T,  t)  dt 


T. 

J 


is  approximated  in  this  fashion.  The  difference  quantities 

rjk  (T  *  V  -  dj(T»' k  =  1’ 2 . ”■ 

for  j  fixed,  are  then  available  for  evaluating  parameters  of  the  random  component 
of  the  impulse  response,  such  as  stationarity,  average  power,  etc. 


i 
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Independent  measurements  of  d(r,  t)  and  D(f,  t)  would  yield  information 
as  to  the  validity  of  the  linear  model,  for  only  in  this  case  are  these  expressions 
related  by  a  Fourier  transform.  Also,  if  corroborated,  the  easier  of  the  two 
measurement  techniques  can  be  selected  and  recommended  for  future  measure¬ 
ment  programs.  (Notice  that  we  can  only  relate  d  and  D,  and  not  h  and  H;  the 
latter  are  noisy,  including  n(t),  and  would  be  impossible  to  relate,  d  and  D  are 
approximately  related  by  a  Fourier  transform  because  they  are  the  means, 
respectively,  of  h  and  H. ) 

The  only  remaining  item  of  importance  in  Figure  1  is  the  additive 
noise  n(t).  This  can  be  observed  with  no  signal  transmitted,  either  determin¬ 
istic  or  random,  its  average  behavior  and  power  determined,  its  low  order 
p.  d.  f's  approximated,  etc.  This  noise  arises  from  water  and  wave  movement, 
fish,  and  other  ships. 

Although  measurement  of  the  transmitted  signal  is  not  directly  related 
to  measurement  of  the  medium  characteristics,  a  few  words  about  it  are  appro¬ 
priate  at  this  point,  in  view  of  the  background  developed  above.  In  order  to 
evaluate  the  amount  of  deterministic  component  s(t)  present  in  the  transmitted 
waveform,  and  least  contaminate  the  results,  and  obtain  the  most  consistency, 
initial  measurements  of  this  quantity  should  probably  be  done  with  a  stationary 
submarine  using  an  absolute  minimum  of  auxiliary  equipment  in  a  quiet  envir- 
onmeiil.  A  long  duration  "close  aboard"  measurement  should  be  conducted, 
repeating  a  basic  signal  (~  100  times).  (The  "close  aboard"  feature  eliminates 
medium  effects;  however  the  measurement  should  be  conducted  far  enough  away 
that  near-field  effects  are  negligible. )  The  total  duration  of  the  run  should  not 
be  so  long,  however,  as  to  allow  any  important  parameters  to  change.  An 
estimate  of  s(t)  is  then  realized  by  averaging  the  100  runs,  properly  aligned 
in  time.  The  difference  between  the  transmitted  waveform  and  the  estimate 
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of  s(t)  it  the  random  component  a(t),  and  can  itself  be  analyzed  to  find  the 
average  power,  degree  of  stationarity,  and  first,  second,  or  higher  order 
p.  d.  f.  's  (probability  density  functions).  Evaluation  of  low-order  p.  d.  f.  's 
may  indicate  what  type  of  process  the  random  component  is,  such  as  Gaussian. 
From  the  above,  among  other  things,  quantitative  estimates  of  the  percentage 
of  deterministic  component  in  the  transmitted  waveform  will  be  available,  and 
decisions  as  to  whether  to  attempt  to  use  the  random  component  in  detection 
and  communication  can  be  made. 

6.  3  Effects  of  Non- Stationarity 

Measurements  of  the  statistical  parameters  of  a  random  process  are 
virtually  always  determined  through  time  averages,  and  the  ergodic  theorem 
employed  to  relate  these  to  corresponding  ensemble  average.  However,  for 
a  non- stationary  process,  these  two  averages  are  not  equal.  It  is  therefore 
of  interest  to  see  what  effect  non- stationarity  has  on  some  of  the  typical  time 
averages  often  computed. 

Consider  first  a  process  in  which  the  non- stationary  character  has  a 
bandwidth  W^;  i.  e. ,  significant  changes  in  statistical  parameters  occur  in  a 
time  l/W^.  Also,  suppose  the  bandwidth  of  the  process  itself  is  W.  Now  if 

JL  «  _L 

W  W  ’ 
n 

many  independent  samples  of  the  process  are  available  in  a  time  constant 
(l/W^)  of  the  non- stationary  fluctuations.  In  this  case,  estimation  of  local 
process  parameters  is  possible  from  a  single  member  function  of  the  ensemble. 
For  example,  an  estimate  of  the  mean,  from  the  j-th  member  function  of  an 
ensemble,  is 
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i  J  x0)  (t)  dt  (T  -  1/Wn) 


2WT 


s-  —  y 

2WT  L 


n=l 


(This  choice  of  T  guarantees  that  x^  (t)  is  quasi- stationary  in  the  averaging 
interval  T, )  If  2WT  >>  1,  (W  >>  W^),  we  have  a  reasonably  good  estimate 
of  the  local  mean. 


If.  on  the  other  hand,  W  and  W  are  comparable,  good  estimates  of 

n 

local  parameters  are  not  possible  from  a  single  member  function.  However, 
good  estimates  of  average  local  parameters  may  or  may  not  be,  depending  on 
the  particular  parameter.  Let  us  illustrate  this  contention.  Consider  a  random 
discrete  voltage,  which  at  time  k  seconds,  can  take  on  voltage  values  k,  k  +  1,  . .  .  , 
k  +  5,  with  equal  probability  (  a  non- stationary  die).  Suppose  we  attempt  to  deter¬ 
mine  the  average  mean*  of  the  voltage  over  an  interval  of  N  seconds  according 
to  the  rule 


N 


k=l 


where  is  a  voltage  sample  at  time  k.  Now 


giving 


*k  =  k+  2 


—  _  N 

y s  t  +  3 


Therefore 


N 


N 


*  s  £  (l‘k-k-f)*  s  £ 


k 


k=  1 


k=l 


The  word  "mean"  indicates  an  ensemble  average;  the  word  "average1 
indicates  a  time  average. 
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But  the  distribution  of  z,  is  uniform  over  the  value*  - 

k  2  2 

and  i*  independent  of  k.  Therefore  a  (z, )  ■  a  ,  and 

K  Z 


5 

2* 


o2(  y) 


1  2 
N  °z 


3 


"  2 


1 

2 


1  3  5 

2  *  2  *  2  * 


Therefore  since  the  variance  of  y  tends  to  zero  as  N  increases,  y  tends  to  the 
average  mean  of  the  random  process  over  any  (large)  interval.  Thus  the  pre¬ 
sence  of  non-stationarity  does  not  preclude  obtaining  some  useful  average  local 
parameters  from  a  single  member  function.  (Of  course,  could  not  be 
evaluated  in  this  manner. ) 


On  the  other  hand,  if  we  attempt  to  estimate  the  average  mean-square 
of  the  voltage  according  to  the  rule 
N 

2 
*k 

k=l 


W  a 


n 


we  find  a  different  result  altogether.  We  have 
x^  =  |  [  k2  +  (k  +  l)2  +. . .  +  (k  +  5)2] 

■  k2  +  5k  +  ^ 


Therefore 


N 


k=l 


N 


w-w  ■  S  I  '  k*  ‘  5k  -  T 1 "  ‘ 


n 


k=i 


where 


_  2  2 
*k,xk*  *k  * 


Therefore 


2  .  .  4  2 

C  V  ’  *k  "  *k 
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But 


=  -g  [  k^  +  (k  +  1  )^  +. . .  +  (k  +  5)^] 

~ 

Combining  this  result  with  the  above  one  for  x^  ,  we  obtain 

2.  .  35  .  2  ,  175  ,  ,  5 

cr  <«kl-  Tk  +  -J-  k  +  79-jJ 


and 


N 


2.  IV  ,35  .  2  175  _  5 

'"'■kZ  <Tk  +  — ^ k+75  36 


k=l 


which  tends  to  infinity  as  N  does.  Thus  the  variance  of  the  estimate  gets 
arbitrarily  large  as  the  averaging  interval  (N  seconds)  increases.  *  Here  we 
would  require  many  member  functions  to  obtain  a  good  estimate  of  the  average 
mean- square  voltage  over  a  long  interval. 


to 


then 


In  general,  when  obtaining  an  average  of  a  random  variable  according 
N 


■  s  1  V 


k=l 


N 


'(y)  =  -7  y 

N2  “ 


a2  (xk) 


k*l 


if 


are  uncorrelated.  Then  if 

1)  a2^)  »  c,  <J2(y) — M 

2)  cr2^)  =  ck,  o2(y)  ~  c 

3)  o2(xk )  =  Ck*,  a  <  1,  o2(y)— ^0 

as  N  increases.  Thus  the  rate  of  increase  of  variance  of  {  x^}  is  of  paramount 
importance  in  the  adequacy  of  the  end  result. 


It  should  be  noted,  however,  that  w/<r(W)  approaches  zero  as  N  increases. 
Thus  the  relative  error  in  estimation  of  ^decreases. 
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In  considering  an  estimate  of  an  average  local  parameter  of  a  process, 
the  merit  of  the  estimate  depends  heavily  on  its  variance.  If  the  variance  is 
small  (in  some  sense),  there  will  be  little  difference  in  estimates  obtained  with 
different  member  functions.  And,  indeed,  this  observation  allows  for  a  relatively 
simple  check  on  the  adequacy  of  a  particular  estimate:  two  estimates  of  the 
desired  parameter  are  calculated  from  two  different  member  functions.  If  the 
estimates  are  reasonably  close  (in  some  sense),  either  estimate  is  adequate. 

If  they  are  widely  different,  the  method  of  obtaining  the  estimate  is  totally 
inadequate. 

We  have  said  nothing  about  stationarity  in  the  above  paragraph.  The 
estimate  will  be  an  approximation  to  a  time  average  of  some  statistical  para¬ 
meter.  This  parameter  may  or  may  not  be  appropriate,  depending  on  the 
problem  under  investigation.  For  example, 


is  an  unbiased  estimate  of  the  average  mean  of  x(t)  over  the  interval  (tj,  t^). 
(This  follows  because  the  ensemble  mean  of  y  is 


2 

where  m(t)  is  the  ensemble  mean  of  x(t). )  If  a  (y)  is  small,  any  two  member 
functions  will  yield  essentially  the  same  value  for  y,  with  probability  near  unity. 
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6.4 


Recommended  Measurements 


In  view  of  the  relative  simplicity  of  the  linear  model  approach  to 
underwater  communications)  it  is  recommended  that  a  concerted  measure¬ 
ment  program  be  conducted  to  verify  the  model  and  ascertain  its  utility  in 
typical  applications.  It  may  turn  out  that  although  the  model  is  correct,  the 
random  fluctuations  in  the  transmitted  waveform,  the  medium,  and  the  en¬ 
vironment  are  so  large  that  satisfactory  communication  and  detection  beyond 
certain  ranges  and  in  certain  geometries  is  impossible,  even  with  the  so-called 
optimum  receiver.  This  fact  of  life,  if  true,  could  be  determined  by  this 
program. 


The  end  products  of  the  measurements  should  be  the  following  items: 


1.  Non-random  signal  component  of  transmitted  waveform  when 
communication  is  taking  place. 

2.  Random  component  of  transmitted  waveform:  mean,  mean 
square,  correlation  function,  second-order  correlation  function,  first, 
second,  and  third  order  p.  d.  f.  's.  (These  quantities  may  vary  with  time. ) 

3.  Non-random  component  of  impulse  response  of  medium. 

4.  Random  component  of  impulse  response  of  medium:  same 
statistics  as  under  2. 

5.  Additive  noise:  same  statistics  as  under  2. 

6.  Non-random  component  of  received  waveform. 

7.  Random  component  of  received  waveform:  same  statistics 
as  under  2. 


8.  Non-random  component  of  medium  transfer  function,  amplitude 
and  phase. 

9.  Crosscorrelation  of  random  radiated  signal  with  received  signal 
(no  deterministic  signal  transmitted). 
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In  order  to  determine  the  quantities  above*  continuous  recordings  of 
the  waveforms  "close  aboard"  and  at  the  receiving  site  should  be  obtained. 
The  pertinent  processing  can  be  done  later.  If  other  statistics  are  deemed 
important  later,  they  can  be  computed  from  the  tape  records. 
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7. 


Discussion 


The  model  of  the  communication  system  suggested  here  is  probably 
the  simplest  one  bearing  any  reasonable  resemblance  to  the  actual  situation. 

How  accurate  it  can  be  made  will  have  to  be  borne  out  by  measurements.  If 
measurements  indicate  that  the  model  is  a  poor  one,  other  possibilities  will 
have  to  be  investigated.  One  simple  alternative  is  a  combination  of  a  linear 
filter  with  memory  in  cascade  with  a  nonlinear  filter  without  memory  (although 
perhaps  time  varying).  However  the  analysis  and  optimization  problems  will 
probably  be  unwieldly  except  by  lengthy  computer  effort.  Very  little  useful 
analytical  work  on  nonlinear  networks  has  been  accomplished  thus  far. 

A  method  of  circumventing  the  whole  model  problem  is  suggested  by 
likelihood  ratio  detection  and  communication.  A  high  order  density  function 
of  a  set  of  samples  is  obtained  for  each  communication  signal  and  for  noise 
alone.  The  decision  as  to  which  signal,  if  any,  was  transmitted  is  reached 
by  computing  the  likelihood  ratios  and  choosing  the  largest  as  corresponding 
to  the  actual  situation.  This  approach  requires  a  great  amount  of  data 
acquisition  and  reduction,  mainly  because  the  number  of  samples  necessary 
to  accurately  estimate  high  order  p.  d.  f.  's  is  enormous.  However  it  benefits 
from  the  fact  that  source,  medium,  and  receiver  need  not  be  separable,  as 
they  were  considered  above,  but  may  be  treated  as  an  entity  with  merely  one 
input  and  one  output. 

A  model  of  a  communication  channel,  similar  to  that  suggested  in 
Section  2  has  recently  been  considered*  and  solved  for  the  optimum  receiver. 
However,  the  relative  amount  of  improvement  afforded  by  the  optimum  receiver, 
as  contrasted  with  the  usual  simpler  matched  filter  for  the  deterministic  com¬ 
ponents,  has  not  been  determined,  due  to  mathematical  intractabilities.  Further 
analytical  work  along  this  line  is  necessary  for  the  underwater  application,  where 
the  random  components  of  signal  and  medium  are  very  substantial  in  comparison 
with  the  deterministic  components. 

*[11,  12] 
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For  the 


The  program  of  study  underway  presently  is  now  outlined. 

*  ! 
linear  time-varying  stochastic  model: 

a)  Relate  the  statistical  properties  of  components  (w^(t)}  of  the 
received  waveform  to  those  of  the  medium  and  transmitted  signal.  These 
properties  will  include  ensemble  mean,  power,  correlation  function,  structure 
function,  higher  order  moments,  and  first  order  probability  density  function. 
Calculate  the  signal -to -noise  ratio  of  the  received  wave.  Investigate  the 
validity  and  effect  of  any  approximations  used. 

b)  On  the  basis  of  the  results  of  a),  decide  which  components  are 
useful  and  what  form  of  receiver  is  desirable  for  detection  and/or  communi¬ 
cation.  This  approach  will  be  limited  to  linear  filters  and  simple  nonlinear 
operations. 

c)  For  a  linear  receiver,  relate  the  statistical  properties  of  the 
output  to  those  of  the  transmitted  signal,  medium,  and  receiver  impulse 
response.  The  properties  to  be  investigated  are  ihe  same  as  listed  under  a). 

d)  Maximize  the  linear  receiver  output  signal-to-noise  ratio  by 

» 

choice  of  impulse  response. 

e)  Evaluate  the  detection  probabilities  at  the  output  of  the  receiver 
for  several  impulse  responses;  choose  those  cases  that  are  most  easily 
realized  first. 

f)  For  known  transmitted  signal  and  receiver  response,  determine 
what  target  characteristics  can  be  deduced  from  the  processed  echo.  Deduce 
what  form  of  sounding  signal  is  optimum. 

g)  For  assumed  statistical  forms  of  medium  and  noise,  deduce  the 
optimum  receiver.  Evaluate  the  performance  of  this  receiver  and/or  approx¬ 
imations  to  it  in  terms  of  signal-to-noise  ratio,  detection  probabilities,  and 
error  probabilities. 
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h)  In  all  the  above,  determine  how  the  receiver  should  adapt  to 
time-varying  changes  in  the  medium  and  noise.  Determine  what  additional 
measurements  are  necessary  to  track  the  changes,  and  whether  the  improve¬ 
ment  is  significant  compared  to  an  average  design. 

i)  Investigate  generalizations  of  the  model  and  results  to  include 
spatial  parameters. 

Almost  all  the  tasks  cited  above  deal  with  the  analysis  of  a  cascade 
of  two  linear  (time-varying)  filters  on  a  transmitted  stochastic  signal.  The 
second  linear  operation  need  not  be  a  strict  match  to  the  transmitted  signal, 
however;  as  such,  it  could  be  called  a  generalized  correlation  receiver. 

That  is,  a  more  general  "multiplier-integrator"  form  is  allowed  than  con¬ 
ventionally.  Results  on  this  configuration  will  constitute  a  very  worthwhile 
addition  to  the  signal  processing  problems  of  several  fields. 
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APPENDIX  A 


M-ARY  WEIGHTED  RADIOMETER 


The  optimum  detection  of  small  (fading)  signals  in  the  presence  of 
noise  has  been  considered  extensively*  in  terms  of  structure  and  signal-to- 
noise  ratio.  Here  we  wish  to  extend  the  analysis  to  communication  and 
accompanying  error  probabilities.  Furthermore,  the  communication  mode 
will  not  be  limited  to  binary,  but  rather,  to  general  multiple  alternative 
(M-ary)  detection. 


To  be  specific, 
over  the  time  interval  (t 
rate  is 


one  of  M  equiprobable  alternatives  is  transmitted 

,  t,  ).  Letting  T  =  t  -  t  ,  the  source  information 
l  b  b  sl 


H  = 


lg  M 
T 


where  lg  x  =  log^x.  The  transmitted  signal,  without  loss  of  generality,  will 
be  denoted  by  Sj(t): 


iu  t 

s^t)  =  Re  {|  L(t)  e  °  } 

where  §  ^(t)  is  the  complex  envelope  of  the  signal,  and  is  non-zero  only  over 
(t^,  t^).  At  the  receiver,  the  "signal"  component  of  the  received  waveform  is 

iu  t 

Re  Uj(t)  G(t)  e  °  } 

where  G(t)  is  the  complex  (low  pass)  medium  (fading)  gain  (see  section  4).  The 
total  received  waveform,  including  noise,  is 

iu  t 

Re  {[£  L (t)  G(t)  +  N(t) ]  e  °} 


The  optimum  receiver  is  an  M-fold  duplication  of  the  basic  weighted 
radiometer,  **  the  only  difference  being  the  pre -multiplication  function  in  each 
branch.  A  block  diagram  of  the  j-th  branch  of  the  M-ary  weighted  radiometer 
is  shown  below. 


*[4] 

** 

Ibid,  p.  49 
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The  output  of  this  branch  is  readily  demonstrated  to  be  given  by 

r.  =  jj[ £  1(t1>G(t1)  +  N(tj)]  [  £  *(V  G*^t2)  +  N*(t2>]  6*<tl)  V  t2>  +*(tl ’V  dtl  *2 

for  all  1  <  j  <  M.  (Infinite  limits  are  allowed  on  the  integrals  since  £  .(t)  is 
identically  zero  for  t  not  in  (t^,  t^). )  $(t)  is  the  complex  low  pass  correlation 
of  the  medium. 


Now  we  shall  make  two  reasonable  assumptions  about  the  transmitted 

signals;  first,  the  bandwidth  of  all  the  product  functions  {£  .(t)  £  (t)}  ,  j  ^  k, 

J  K 

will  be  much  wider  than  B,  the  "spread"  of  the  medium.  Since  B  =  1/t  ,  where 
tq  is  the  coherence  interval  of  the  medium  fading,  this  first  assumption  is 
equivalent  to  many  "wiggles"  of  the  product  functions  in  time  tq.  Secondly, 
the  transmitted  signals  will  have  flat  amplitude  distribution  in  time;  i.  e. , 

|  £  j(t)  |  =  '  4a  <  t  <  V  aU  j' 


where  S  is  the  average  power  of  the  received  signal, 
a 


Since 


sa=  \  IfijWCHt)!2  =  J  |£j(t)  I2  |G(t)  |  2 


we  must  have 


|  G(t)  |  2  =  1 

which  is  no  more  than  a  choice  of  normalization.  More  generally,  we  have 
♦<tl-t2)  = 


*Ibid,  pp.  50-54. 
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the  correlation  function  of  the  medium  fading.  Using  the  above  assumptions, 
we  have 

s, 

rl  s  II  [2S&G(t1)  +  e*1(t1)N(t1)]  [2SaG*(t2)  +  ei(t2)  ♦Vj  -  t2)  dtj  dt2 

t 

a 

and 

•b 

r.  S'  jj  Hit,)  N*(.2)  £*(»,)  S.(«2)  ♦*(.,  -  «2),  j  >2. 

t 

a 

Utilizing  the  Gaussian  and  white  properties*  of  N(t),  we  find 

77  =  (2  S  )2  I  +  8  S  TN 
1  a  2  ad 

where 

T/2  T 

-  JJ  I  ^*1  "  V  I  2  dtl  dt2  =  J  (T  *  I  u  I  >  I  I  2  du 

-T/2  -T 

and  N.  is  the  double-sided  noise  power  density  level  (N  =  2N  ). 
d  -  o  d 

Also, 

77  =  8  S  TN  ,  j  >  2. 
j  a  d  — 

Continuing,  it  is  tedious,  but  not  difficult,  to  show  that 

j-  i 

o2(r  )  =  (2  S1)2  (4  Nd!2  I2,  j  >  2 
'  «2(rl)  =  (2Sa)2(4Nd)2  I2  „ 

*[3] 
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where 


SI,  /  S  \  I. 

“  Nd  .  *2  [2Nd  )  Z2 

T/2 

“3  ■  jlf  +(tl -*2>  *(,2  -*3>  ♦VV  d,l  dt2  dt3 

-T/2 

T/2 

J4  =  "t2)  *(t2  "t3)  *(t3  “  V  +(t4"tl)  dtl  dt2  dt3  dt4* 

-T/2 

Now  if  T  >>  T  ,  many  independent  segments  of  input  signal  will  be  added  together 
(with  noise)  to°form  the  variables  {r.}.  Therefore  {r.}  will  satisfy  a  joint 
Gaussian  distribution.  The  probability  of  correct  detection  is  then  readily  demon¬ 
strated  to  be* 

p  =  f  e*E(-xZ2,i  [ay  +  g]  dy 
C  J  J*i 


where 


N 


4. 

J  (T  -  |  u|  )  |  <)>(u)  |  2  du 


-T 


1/2 


Now  since  the  Gaussian  approximation  holds  only  when  T  >>T0»  we 


S 

IT 

O 


oo 


J  T  |4>(u) 


|  2  du 


-80 


J 


1/2 


For  a  linear  correlation 


4>(u)  =  1  -  .  |  u|  <  2T£ 

o 

fA  S  T  S  t 

3?i  — a—  -J—2 

H  V3  N  N 


'[9] 
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Also  a  ~  1  if  S  t  /N  <<  1.  Therefore, 
a  o  o 


S  T 
a  o 


N 

o 


dy. 


This  result  is  used  extensively  for  plots  given  in  the  main  body  of  the  report. 
For  M  =  2,  a  special  case  results: 


If  the  source  data  rate  is  kept  fixed,  and  M  increases  without  bound, 
we  can  use  the  result  that 


lim  I  x  +  y  yjl  In  M'  — >1,  all  x,  if  Y  >  1 

M— $oo  l 


to  show  that 


S 

P  — $  1  as  M— )  oo  if  H  <  lge  rr— • 
c  N 


Thus,  analogous  to  the  non-fading  case,  zero  error  probability  can  be  obtained 
in  the  limit,  with  a  non-zero  source  rate.  Thus,  the  capacity  of  a  fading 
medium  (infinite  bandwidth)  is  bounded  by 


This  relation,  and  the  others  in  this  appendix,  hold  if 
S  80 

rf-  J  2  du  <<  l. 


Furthermore,  in  this  case,  and  an  S  =  S  power  limitation,  this  is  the  optimum 
obtainable  performance.  ^ 
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APPENDIX  B 

MULTIPLE  ALTERNATIVE  RECEIVER 


This  appendix  will  be  concerned  with  the  derivation  of  the  error 
probabilities  of  a  multiple  alternative  receiver,  where  the  same  alternative 
(message)  is  repeated  sequentially  for  the  entire  message  interval  T.  The 
derivations  are  conveniently  broken  into  two  main  parts:  the  first  part  deals 
with  the  case  of  no  peak  power  limitation,  whereas  the  second  part  considers 
a  peak* equal-average  power  limitation. 


No  Peak  Power  Limitation  (S  >  10  S  ) 


Without  a  peak  power  limitation,  it  might  appear  at  first  that  the 
signal  energy  should  be  packed  into  a  time  interval  short  compared  with  the 
coherence  interval,  t  ,  so  as  to  avoid  noncoherent  fading  throughout  the 
message  interval  T.  However,  an  immediate  and  obvious  drawback  to  this 
proposed  technique  is  that  communication  is  too  heavily  dependent  on  the 
instantaneous  fading  of  the  medium.  A  much  more  desirable  alternative 
would  be  to  split  the  total  signal  energy  into  several  pulses,  each  of  duration 
less  than  tq,  and  then  add  the  components  after  detection.*  This  method  per¬ 
forms  some  averaging  over  several  fading  intervals.  Specifically,  a  possible 
partitioning  of  the  signal  energy  is  depicted  in  the  figure  below. 


There  are  n  pulses  in  the  message  interval  T,  of  equal  height  A.  For  an 
average  power  constraint  S^,  we  must  have 


S 

a 


le.w  G(t> 

j 


I  dt 

t 

a 


nT 


Therefore, 


A 


2 


2  S 

a 


A 

T 


Here  we  are  allowed  to  vary  A  and  7  as  we  please,  and  choose  A  so  as  to 
satisfy  the  power  constraint.  The  received  waveform  is 

e ‘“s'  } 

We  define  R.,  as  the  envelope  of  the  output  of  the  filter  matched  to  the  j-th 
message  during  the  k-th  pulse.  Here,  1  .£  j  <  M,  l<k£n.  Then 


Re  {[S^t)  G(t)  +  N(t)  ] 


{k-l)A+T 

j  e.*(t)  [ e  !<t)  G(t>  +  N(t>]  dt 

(k-l)A 


all  j. 


Now  for  constant  amplitude  orthogonal  signals  over  the  interval  of  duration  r, 


(k-l)A+r 

(k-l)A+r 

Rlk  = 

2  S  — 
a  t 

^  G(t)  dt  + 

j  £  L*(t)  N(t)  dt 

(k-l)A 

(k-l)A 

T 

T 

s 

2  S  — 
a  r 

J  G[  t  +  (k-l)A  ] 

dt  +  j  e  j*(t)  n  [t  +  (k 

o 


o 
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where  we  have  set 


i i  [t  +  (k-l)  A]  =  €  j(t).  0  <  t  <  t  ,  all  j, 
corresponding  to  the  repeated  signal  assumption.  Also 
T  T 

j  i.*(t )  €  x(t)  G[t  +  (k-l)A]  dt  +  jg.*(t)N[t+(k-l)A]dt 

o  o 

Now  we  consider  adding  the  v-th  powers  of  the  envelopes  of  the  j-th  message, 
to  form  the  decision  variables: 


k=l 


The  source  information  rate  is  given  by 

H=  itM  «  li_M 

T  nA 


For  fixed  A  and  H,  as  M — *eo,  so  also  must  n.  Therefore,  each  of  the  yj 
variables  approaches  a  Gaussian  variable,  and  we  can  compute  the  probability 
of  error  providing  we  can  evaluate  y.  and  o^(y.).  We  start  with 


But 


n 


k=  1 


where 

T 

zjk=  {  6  j(t)  G  [t  +  (k-l)A  ]  +  N[t  +(k-l)A]  }  dt,  all  j 

o 
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Due  to  the  linear  operation  on  G  and  N,  both  of  which  are  Gaussian,  x  and 

y.,  are  Gaussian.  Furthermore,  J 

7jk 


Z..  Z#  =  0 

jk  /m 

and,  after  simplification, 


ZjkZ*/n,  *  H  s,*(tl)  W  £  dt; 


+  4Nd(A  T»  WjJ 


using  the  orthogonality  of  the  signals.  In  order  to  make  further  substantial 

progress,  we  assume  that  A  >  4tq  and  T  ,5  J  To'  Pky»ically  this  corresponds 

to  uncorrelated  fading  between  pulses  and  negligible  fading  within  a  pulse, 
respectively.  With  these  assumptions,  there  readily  follows 


Z.,  Z*. 
jk  im 


..2  .2 
=  (A  r) 


6km  6jl  6l  1 


(A  T)  4Nd  ‘km 


This  formula  indicates  that  all  the  statistical  variables  in  different  pulses  are 
independent  i  Therefore  we  need  confine  our  analysis  to  a  single  pulse  and 
later  scale  the  means  and  variances  by  n.  Dropping  (unnecessary)  second 
subscript,  we  have,  expressly, 


Z,  Z*  *  (A2t)2  +  (A2t)  4N 
li  a 

=  o,  j  H 

z.  z*  =  (A2T)  4N  ,  j  ji  1. 

J  J 


Therefore 


xi  =  y. 


/A2  x2 

(A  t) 


+  (A  t)  2N  • 
a 
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2  2 
x .  -  y  . 
3  7  3 


(A2t)  2Nd  *  a2  ,  j  i  2 


y.  =  x.y.  =  0 
J  3  3 


and 


2  - 1 

p(x2»  Yj)  =  (2 to  j  )  exp 


2-1 

P^Xj ,  y^)  =  (2ira2)  exp 


X1  +  Vl 


2  a 


1 


2  ,  2 
x  .  +  y  . 

_J _ _J_ 

2  a? 


,  j>  2. 


Therefore  since 


R.V  =  (x2  +  y2) 
J  3  7  3 


v/2 


v/2 


Rj  3  II  Uj  +yj}  (2waj)’1  exp  - 


-  oo 


2  2 

Vyi 

2  a2 
J 


dx.  dy. 
J  J 


■2“/2  •;  -  &  ♦ 1 


Also,  immediately, 


RjV=  2VffjVr<v+1> 


Therefore 


yj =  n2''/2  ,’vj  r  (f  ♦  ■] 


and 


o  (y.)  =  n 

J 


R2v  -  R v 
J 


;] 
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„v  2 
»n 2  a 


T  [r-i-r2  |l+l)] 


Therefore  the  signal-to-noise  ratio,  upon  which  error  probabilities  critically 
depend,  is 


’± 


_  2  fv 

n  r  |f  +  1 

r(v  +  i)  -  r2  j  +  i 


and  the  corresponding  probability  of  correct  decision  is 


f  exp  (-x2/2) 

_  M- 1 
w 

1 

-i 

_a(y2) 

a(y2)  J 

Simplifying,  we  have 

<r(y2)  ~ 


'/2 


,1  + 


.2  1 

A  T 

4N , 


a  (  l  + 


V  \ 

No  I 


>/z 


dx 


and 


*1  ~  y2 
o-(y2) 


= 


or  =  3 

V 


where 


a  * 
v 


r  2  +1 


~\Jr(v  +i)  -  r2  ||+  ij’ 

In  order  to  obtain  numerical  results,  let  us  now  choose  A  a  l/H;  that  is,  let 
the  pulse  repetition  period  equal  the  bit  interval.  Then  since 

„  lg  M  lg  e  In  M 
T  nA 


we  must  have 


n  *  lg  e  In  M 
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! 


Then 


v/i 

8  ®  Ylg"e  |l  +  |  -1  av  -y/iiT 


exp  (-x  /2)  _  M- 1 


I  •  i\v/2  1 

(■^1  H 


This  is  the  general  expression  for  the  error  probability  as  a  function  of  M, 
signal-to-noise  ratio,  v ,  and  information  rate.  We  now  use  the  fact  that 


I™  ®M  [y  +  y  V 2  In  M  ]  =  1.  all  y,  if  y  >  1, 


to  show  that 


p  — ■»  1  as  M  - r«o  if 

c 


S  /N 
a  o 


The  right-hand  side  of  this  last  equation  is  a  lower  bound  on  the  capacity  of  a 
fading  channel.  The  representative  numbers  given  in  the  table  below  indicate 
that 


1  1.5  |  2  2.5 


4  6  I  8  I  10  I  12  I  oo 


S  /N 
a  o 


.  621.  75  1.  85  1.93  |  1.  0  |  1.  1  |  1.  2  1 1. 23  1 1. 23  1.  22  1.  0 


the  capacity  of  a  fading  channel  is  approximately  lg  e  =  1.  443  ^/N^, 

just  as  though  no  fading  were  present.  This  topic  deserves  much  further 
analysis,  particularly  regarding  the  applicability  of  the  Gaussian  form  as  n  and 
M  simultaneously  tend  to  infinity.  However,  thfls  above  results  hint  that  proper 
use  of  a  fading  channel  can  lead  to  almost  the  same  performance  as  in  the 
absence  of  fading. 
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) 


B.  2  Peak-Equal-Average  Power  Limitation  (S  ■  S  ) 

_ _ _ E _ £_ 

2  Uaing  the  same  notation  as  before,  in  this  case  we  must  have  A  *  r, 

and  A  =  2S  ;  therefore, 
a 


z..  z*~  = 

jk 


$$  V'i>  W £  i(ti>  £I  'V  ♦  t’r  V  )  d'i  d,2 


+  4N  (2S  t)  6  6.  . 

d  a  km  )£ 


Now,  restricting  t  to  be  smaller  than  —  r  ,  we  get 

Z  o 


zjk  Z*J| m *  <2V>  sji^i+«k-“w+4Nd(2V>skmsie 


I 


Now  we  introduce  the  restriction  v  -  2.  In  this  case, 


and 


n 


*  I  V2  ■  I  V  z] 


* 

jk 


k=l 


k=l 


Also, 


V  t  zikZ*k=  ",zv»  1(2  V1  V  +  4V 


k=l 


n  n 


z  z*  z  z* 

jk  jk  jm  jm 


k=l  m=l 


i. 


n  n 

H 

k=  1  m=l 


z.,  z*  z.  z.  +  z..  zf  z*  z. 

jk  jk  jm  jm  jk  jm  jk  jm 
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(4 


since  Z  .  is  Gaussian  with  zero  mean.  Using  the  expression  for  Z..  Z*.  above 
jk  j*  4m 

we  get  Z..  Z*  Z*  Z.  -  (2  Sr)2  [2  Sr  6 +4N  ]2 
•  jk  jm  jk  jm  a  a  jl  d 

and 


zTz T~  Z*  Z.  =  (2  S  t)  [2  S  r  +  ((k-m)r)  +4N.6.  1 

jk  jm  jk  jm  a  1  a  jl  d  kmJ 


Thus, 


r2  *  (2  S  r)2  j  n2  [2  S  r  6..  +  4  N,  ] 
j  a  C  a  jl  d 


n  n 

+ 

k=  1  m=  1 


I2  V*  «k-”>T>sji+4Nd8to„)2j 


(2Sar)2  [„2  [2S>T«.,t4Nd]: 


n  n 

+  (2SaT)2  ^  ^  <t>2((k'm)T)  6jl 

k=  1  m=l 

+  1 6  n  S  r  N ,  6. , 
a  d  jl 

+  16  n  j 

Now,  making  use  of  the  approximation, 

nT 

[  (nr-x)  f(x)  dx. 
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we  get 


(2 


S  T) 
a 


l 2  V  !ji 


+  4Ndl 


nr 


2<2Sa»2  5 


(nr-x)  $  (x)  dx  •  6 


jl 


+  16  n  S  t  N,  6.,  + 
a  d  jl 


16nN2| 


.2  2  — 2 
”Q\  "  VTJ 


*  n(2  S  r  N.)  •  16 
a  a 


nr 


.  ,-jf  -fc  5  <»T-xM>2(x)  dx 

LV  al 


S  r 
a 

N. 


6. .  +  1 

Jl 


Thua, 


This  expreasion  is  maximized  by  letting  r  *  T  /2  (the  largest  value  allowed  by  the 

S  T 
a 

analysis)  thus  giving  (for  <<  1), 

o 


exp 


„M-1 

v 


S  T 
a  o 

2N 


o 


dx 
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APPENDIX  C 


DELAY  CORRELATION  SYSTEM 


The  blotfk  diagram  of  the  system  to  be  considered  here  is  depicted 

below. 


In  this  system,  the  quantity  D  is  compared  with  a  zero  threshold: 

D  >  °  }  tran,n,i'"d- 

1 

To  evaluate  the  performance  of  this  system,  we  will  derive  R  =  — j-., 

ffD 

which  will  provide  the- bit  error  probability  for  Hr  <<  1.  In  terms  of  the 
received  signal,  v(t),  and  the  linear  filter  impulse  response,  h(t),  the  decision 
variable,  D,  is  given  by 


=  y  y<t)  dt 


where  y(t)  =  x(t)  x(t-T) 


and  x{t)  =  ^  x^(t) 

j=l 

t 


yi»  ■  v  I 


v(u)  h{t-u)  du 


(j-l)T 
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where  it  ie  aiiumed  that  h(t)  =  0  for  t  <  0  and  t  >  r. 

and  Fj(t)  =  1  jr  <  t  <  (j+i)r 

-  0  otherwise 

Substituting  the  explicit  form  of  y(t)  in  the  first  equation,  we  obtain 
n-1  (j+l)T  t  t 

-Z  I  dt  y  dUj  ^  du^  ''(u^)  hft-Uj)  hft-u^) 

j=l  jr  jr  jT 

n-1  (j+l)r 

=  2  I  *  A,W  BjW 

j=l  jT 

t 

where  A^(t)  =  J  R-e  { or t) }  =  ^  du  v(u)  h(t-u) 

jT 

t 

Bj(t)  =  j  Re  {«t)}  =  §  du  v(u-t)  h(t-u) 

jT 

i2irf  t  p 

and  ar^t)  =  e  °  J  V(u)  H(t-u)  du 

jT 

i2irf  t  p 

0.(t)  =  e  °  J  V(u-t)  H(t-u)  du, 
jT 

and  V(t),  H(t)  are  the  complex  envelopes  of  the  received  signal  v(t),  and  filter 
impulse  response,  h(t),  respectively. 
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Now  for  the  mean  value  of  D, 
n-1  (j+l)T 

5=  £  j*  dt  A  (t)  B  (t) 

j=l  jT 


But 


Now, 


A.(t)  B  (t)  =  i  Re  (a.(t)  0.(t)  +  'o.(t)  CTt)T 
J  J  8  J  J  J  J 


*1  *2 


_  i2irf  (t.  +  t  )  p  «  _ 

ajl(tl)  0j2(t2)  =  e  °  J  dui  J  V(Ul)V(u2-r)  IKtj-UjJH^-Uj) 


v  hr 


But,  since  V(tj)  V(t2)  =  0  for  all  tp  t2  when  V(t)  is  Gaussian, 

VV  W  - 0 

Consider  next  <*jj(tj)  &p*^2^  *2’  ^1’  and  $2  8enera^»  since  we  will  need 

2* 

the  general  result  later  for  D  ): 

“•Wvl  c  r 

\  J  dul  J  du2  *(urU2+T)  ^(Ul)  ^*(U2_T) 


WW  =  e 


ijr  v 


HUj-Uj)  H*(t2-u2) 


J.  A. 

+  UNo)  I  dU!  J  du2  6(u1-u2+T)H(tl'Ul)H*(t2'U2) 


JiT  j,T 


I 


See  [2],  p.  53. 
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where  $(u)  ■  G(t)  G*(t-u).  Since  N(t)  ie  derived  from  white  Gaussian  noise.  * 


N(tj)  N*(t2)  =  2  Nq  6{t j  -  t2). 
and  N^  is  the  single  sided  noise  power  density. 


Now,  setting  jj  =  j2  and  tj  ■  t2,  we  obtain 
t  t 

«j(t)  &*{t)  =  J  dUj  J  du2  ♦(uru2  +  T)  *(u2-t)  Httj-Uj)  H*(t2-u2) 

jr  jr 

We  now  let  (for  the  remainder  of  the  appendix), 

|£<t)  j  =  yJzT  0  <  t  <  T 

=  0  otherwise 


and 


H(t)  a  F^t) 


so  that 


n-1  (j+l)r 

r-s  r\ 

t 

t 

D=(2Sa)  2  J  * 

s  iau 

1  J  du2  R*  (♦<uj"u2+t))  F^t-Uj)  F^t-u^ 

j  =  l  jT 

jT 

jT 

/,c  \  n-l(j+l)T 

t 

t 

idui 

J  du2  Re  {(Huj  -  u2  +  t)} 

j=l  jr 

jT 

jT 

Now,  using  the  identity 
a+b 

^  f(x-y+c)  g(x-y+d)  dxdy 


=  J  (b  -  |  u  |  )  f(u+c)  g(u+d)  du. 


-b 


Ibid,  p.  55. 
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we  get 


t  t  t-jT 

j  dUj  ^du2  +(u1-u2+t)»  ^  du<Mu+r)  [t  -  jr  -  |u|]  . 
j  t  jr  -(t-jx) 


Substituting  this  integral  in  the  expression  for  D,  and  transforming  the 
variable  of  integration,  t,  each  term  in  the  sum  can  be  made  independent  of  the 
index  j,  giving 


T  X 


D  =  i  (n-1)  S&  ^  dx  ^  du(x-u)  Re  {♦(r+u)  +  ♦  (r-u)  } 
o  o 

The  region  of  integration  is  the  u,  x-plane  is  shown  in  the  diagram. 


T  T 


Integration  over  this  region  can  be  performed  in  the  form  \  du  \  dx;  using 
this  order,  D  becomes  **  J 


o  u 


T  T 


D  =  i  (n-1)  ^  du^dx  (x-u)  Re{  $  (r+u)  +  ♦  (r-u)} 


o  u 


r-u 


*  ^  (n-1)  S&  j  duRe  {  $  (r+u)  +  $  (r-u)}  ^  dy  •  y 
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■  —  (n-1)  S&  ^  du(r-u)^  Re  {  $  (r+u)  +  $(t-u)} 
o 

=  i(n-l)  Sa  «Wr) 

T 

where  4*  (r)  =»  j  du  Re  {  ^  (2r-u)  +  $  (u)} 
o 

For  the  second  moment  we  have 

_  n^l  n-1  °1+1)T  V)T 

°z  - 1  l  1  d,i  hwwvw 

il-1  j2-l  jjT  j2T 

But 

Ww??v? 


=  i  Re  {a.^tj)}  iRef^.^tj)}  \  Re  {or j2(t2)>  j  Re 


s  ~T  Re  ^1S1  +al  81  *>  Re  <a202  +a2B2  *} 

2 

=  Re  («  jSj  *  a2B2*  +  a  lBl  a  2*  B2*  +  a  1B1*  a2  *  B2 ^ 
2 

where  a ^  ■  a  ( t^),  and  0^  *  Bj/}t^  ~  l,Z 

(Averages  of  all  terms  not  containing  an  equal  number  of  conjugated  and 
unconjugated  factors  are  zero,  since  the  ar's  and  S's  involve  the  same 
Gaussian  process,  V(t)  ). 
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1 


Now, 


a .  8 ,  *  a  «.  S-  B  a  ,  8,*  <*  _  8»*  +af,8.+  a.  8  * 

1  pi  2m2  1^1  2  2  1  2  21 

a  8  a  *  S  *  =  a  <y  *  8  8  *  + «  8  *  a  *  8 

1  P1  2  p2  12  °1  °2  1  p2  2  1 

a  8  *  a  *  8  =  "a  a  *  8*8  +  a  8  *  a  *  8 

1  P1  2  2  12  1  2  11  2  2 


In  the  remainder  of  this  appendix,  we  shall  assume  that  the  signal-to- 
noise  ratio  at  the  output  of  the  linear  filter  is  small,  so  that  only  terms 
involving  are  significant.  With  this  assumption,  we  have,  from  the 
expression  derived  earlier, 


Similarly, 


*1  *2 

iZ’rf  Jt.-  t  )  f  r 

(2No)  ]  dUl  j  du26(Ul-u2  + 

V  j2T 


al“2* 


»  *2* 


i27rf0(tr  t2) 


*1  ‘2 


j  dUj  j  du2  4>(uru2) 


V  J2T 


*1  *2 


+  2  N 


j  dui  j  du2  6(uru2) 


v  hT 


(2  Nq)  e 


i2lrfo(tl-  V 


S  ‘2 


j  dul  \  du2  6  (U1-  V 


J1T  >2T 


I 


I 


Therefore, 


*1  Vtt2  s2*  ■  (2  No)2  iljj  6(V  wl  +  t)  6(u2“  w2  +  T) 

+  (2  No)2  iiH  6(V  U2  +  T)  6(w2"  W1  +  T) 

*1  *2 

Whe”  1151  *  il  dUl  dWl  ii  dU2  dW2 


>l7 


V 


[ 


But  both  of  these  terms  are  zero,  since  |  u^  -  w^  |  <  t  for  £-  1,2.  In 
the  first  term,  the  6  -functions  are  zero  unless  |u^  -  w^  |  >  T,  and  in  the 
second  term,  at  least  one  of  the  6  -functions  is  zero  unless  |  u^ -  |  = 

|  u^-  t  -  w^-  T  |  *  |  u^-  w^-  2r  | ,  in  which  case  |  u^  -  w^j^  r. 


I 


Similarly, 


«1  0i  a2*^2* 


*  ®  1  a2*  S1S2*  +  “1*2*  Q  2*  81 


r 

! 


=  (or  j  or  2*)  +0 

i2irf  2(t,  -  t,) 

=  (2  N  )2  e  ° 
o 


2  iiil  S<W1'  W2> 


i2rrf  2(t  -  t  )  ? 

s(2Nc>  e  (t2"jT)  ifjl=j2=j 


t 

I 


and 


«  1  81*  “2  S2 


=  0 


otherwise 


a  1  °  2*  81*  02  +  “l8l*  a  2*  02 


I®  i«2*  '  +  ° 

(2  No)2  V  jT)2  if  jl  *  j2  “  j 

0  otherwise 


-77- 


Therefore. 


AJ1**1*  BJ1**1*  AJ2**2*  BJ2**2*  ’  7  R"  1 


(2  n/  (t2-  jr>2 


1  +e 


i2irf  *  2(t,  - 
o  1 


■■’]} 


otherwise 


Substituting  in  the  expression  for  D  ,  we  get 
_  n-1  (j+l)r  (j+l)r 

°2  r  1  5  dtl  5  dt2  7"  (2  V2  '*2-  I  1  +  “■  4rio«'l  -  *2> ) 

j=l  jT  jT 

ui  *  -  JT'  ^s1'  2  =£> 

2  n- 1  r  t 


=  I  j  dUl  j  dU2  U2  [  1  +  008  4,rfo(ul  •  u2)  ] 

j=l  O  o 


N 


o 

.5 


(n-1)  —  +  \  cos  4irf  u,  du,  1  u2  cos  4irf  u_du_ 

3  J  olljc  oZZ 


N  2  (n-1)  t4 
o 

3  •  2 5 


j 


+  \  sin  4*01  u.du. 

oil 


j": 


,  sin  4irf  u0du 
2  o  2  2 


if  f  is  much  larger  than  —  . 
o  r 


Now,  combining  the  results  for  D  and  D  ,  using  the  latter  as  an  estimate  for 
2 

oD  (valid  for  low  signal-to-noise  ratio),  we  obtain 


.5(1 


2  „  2  ,2, 


_  D  ~ 
RS  “  = 


3*  2  (n-1)  Sa  *  (r) 

Nq2  (n-1)  t4 
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3  (n-1)  S  2  *2(t) 
a 

"  -  M  2  4 
2  N  t 

o 

T 

Substituting  n  ■  -  .  this  can  be  rewritten 

T 


R?  3  (1.r)±!M. 

t  r 
o 


S  T 
a 

N 


S  T  1 
a  o 

N  ■ 


and  if  T  >>  r,  then 


S  T 


-  3  iT(t)  _a 


2  5 

T  T 
O 


N 


S  T  1 

_a _ o 

N 


The  triangular  form  of  $  (u)  will  now  be  examined;  i  •  *  *et  <j>  (u)  =  1  -  1 2T 

Then 


u 


•Mr)  =  4  t3  1  - 


2t 


R*  f  f  i- 


,2  /ST 


2t 


S  T 
a  o 

N 


2t 

o 

This  expression  can  be  maximized  by  setting  r  -  - 
we  get 


Making  this  substitution 


~  2 


opt. 


fs  T\  is  t 
a  o 

N 


N 


Now,  since  this  is  a  binary  system  with  S  =  Sa>  we  have  T  *  ^  ,  and  therefore 
the  quantity  R  may  be  written  ^ 


r-  II!4 


s  i  1 

Sa  H 


N 


S  r  ' 
a  o 

N 
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For  (Ht  )  >>  1,  the  decision  variable,  D,  is  Gaussian,  and  the  bit  error 
o 

probability  may  be  written 


=  *  (-  /R1). 

For  t  =  2/3  r  , 


p2  =  .  (-<2/3, 2  P  ynr] 


where*  as  before* 
x 

•w  *  J 


u 

2 


du. 
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